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ABSTRACT 


\ 

^ A maximum  likelihood  estimation  program  was  applied  to 
flight  data  for  Princeton's  Variable-Response  Research 
Aircraft  to  determine  its  primary  stability  and  control 

t 

derivatives.  The  control  derivatives-  for  the  side-force 
surfaces  and  the  rudder  were  of  special  interest.  The 
effects  of  measurement  noise  and  process  noise  on  parameter 
identification  also  were  studied. 

This  investigation  shoved  that  the  maximum  likelihood 
estimation  program  used  identifies  derivatives  vhich  produce 


close  fits 

of  the 

■easured  time  histories.  Standard 

deviation 

£ of 

the  derivatives  computed  from  several 

time 

histories 

indicate  the 

quality  of  the 

estimates. 

The 

reduction 

in 

standard 

deviations  vhen 

estimates 

vere 

separated 

by 

type  and 

direction  of  control  input 

ti  me 

history 

used 

indicates 

that  derivative 

estimates 

are 

affected 

ty  the  assumptions  inherent  in  the 

analytical 

model 

and  the  signal 

-to-noise 

ratios  of  the  data. 

The  method 

used 

for  identifying  highly  correlated  derivatives  also  affected 
the  estimates  obtained.  The  final  set  of  derivatives 
determined  in  this  research  produced  a good  fit  of  the 
measured  data  and  several  of  the  derivatives  agreed  veil 
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ABSTRACT 

A maximum  likelihood  estimation  program  was  applied  to 
flight  data  for  Princeton's  Variable-Response  Research 
Aircraft  to  deternine  its  primary  stability  and  control 
derivatives.  The  control  derivatives*  for  the  side-force 
surfaces  and  the  rudder  were  of  special  interest.  The 
effects  of  measurement  noise  and  process  noise  on  parameter 
identification  also  were  studied. 

This  investigation  showed  that  the  maximum  likelihood 
estimation  program  used  identifies  derivatives  which  produce 
close  fits  of  the  xeasured  time  histories.  Standard 


deviation 

s of 

the  derivatives  computed  from  several 

time 

histories 

indicate  the 

quality  of  the 

estimates. 

The 

reduction 

in 

standard 

deviations  when 

estimates 

were 

separated 

by 

type  and 

direction  of  control  input 

ti  me 

history 

used 

indicates 

that  derivative 

estimates 

are 

affected  ty  the  assumptions  inherent  in  the  analytical  model 
and  the  signal-to-noise  ratios  of  the  data.  The  method  used 
for  identifying  highly  correlated  derivatives  also  affected 
the  estimates  obtained.  The  final  set  of  derivatives 
determined  in  this  research  produced  a good  fit  of  the 
measured  data  and  several  of  the  derivatives  agreed  well 
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with  analog  matching  derivative  estimates. 
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1.  INTRODUCTION 


Identification  of  aircraft  parameters  is  becoming 
increasingly  important  in  applying  electronics  technology  to 
aircraft  design  for  a number  of  reasons.  Modern  control 
design  techniques  rely  on  accurate  aircraft  aerodynamic 
models  for  best  results.  computer  analysis  of  aircraft 
stability,  handling  qualities,  and  mission  effectiveness 
also  depends  on  accurate  aerodynamic  models.  Parameter 
identification  from  flight  test  data  offers  the  advantage 
over  analytical  and  wind  tunnel  results  of  determining 
stability  and  control  derivatives  of  the  actual  aircraft  in 
the  flight  regime  of  interest.  The  mechanization  of 
parameter  identification  and  the  application  of  this 
procedure  to  Princeton’s  Varianle-Fesponse  Research  Aircraft 
(VRA)  is  the  subject  of  this;  thesis. 

The  VRA  is  a specially  modified  Ryan  Navion  aircraft. 
Several  studies  have  been  made  to  determine  its  stability 
and  control  derivatives  with  varying  success,  including  mind 
tunnel  tests  of  the  full  scale  aircraft  (Sef  1),  application 
of  analytical  methods,  analog  matching  of  test  flight  data 
(Ref  2)  , and  maximum  likelihood  estimation  from  test  flight 
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data  (Ref  3)  . This  study  will  make  use  of  the  Bodified 
Maximum  Likelihood  Estimation  program  developed  by  Haine  and 
Iliff  (Ref  4)  and  similar  in  concept  to  the  routine  used  in 
Pef  3. 

1.1  Thesis  Summary 

The  objectives  of  this  thesis  include  the  following: 

. Identification:  of  the  primary  longitudinal  and 
lateral-directional  stability  and  control  derivatives 
for  the  VRA  in  its  current  configuration  at  a cruise 
condition  of  105  KIAS.  . 

. Identification  of  control  derivatives  for  the 
modified  rudder  and  the  irecently  added  side-force 
surfaces. 

. Estimation  of  the  effects  of  measurement  and 
process  noise  on  parameter  identification. 

In  addition,  the  results  obtained  will  reflect  the 
ef f iciency^  and  accuracy  of  the  computer  program  used. 

A total  of  28  control  input  time  histories  from  2 test 
flights  were  used.  The  twenty  lateral—  directional  time 
histories  included  aileron  pulses,  rudder  doublets,  and 
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side-force  surface  pulses  and  doublets.  Estimation  of 
longitudinal  derivatives  was  restricted  by  the  availability 
of  only  8 suitable  elevator  doublet  time  histories.  Dynamic 
measurement  errors  in  the  angle  of  sideslip  and  velocity 
instruments  are  documented;  poor  lateral  accleration 
measurements  and  possibly  erroneous  angle  of  attack 
measurements  also  are  noted. 

The  need  for  an  improved  analytical  model  is 
demonstrated  by  the  seeming  ly. nonlinear  responses  in  roll, 
pitch,  and  yaw  to  control  inputs.  The  large  control 'inputs 
required  by  signal- to- noise  ratio  considerations  add  to  this 
nonlinear  aircraft  behavior. 

Improved  control  input  design  is  needed  to  take  full 
advantage.^  of  the  maximum  likelihood  estimation  tecnnique. 
Control  inputs  that  failed  to  excite  all  'modes  of  motion 
resulted  in  large  variances  in  the  estimates  of  derivatives 
related  to  these  modes.  This  result  is  demonstrated  by  the 
poor  time  history  matches  for  aileron  inputs.  The  VBA  now 
has  the  capability  for  precise  control  inputs  through  its 
microprocessor  Digital  Plight  Control  System.  Tailored 
inputs  would  make  the  best  use  of  existing  controls  in 
exciting  all  modes  of  interest. 
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1 . 2 Organization  of  Thesis 

The  remainder  of  this  thesis  is  made  up  of  5 chapters 
and  3 appendices.  Chapter  2 discusses  the  theory  of  maximum 
likelihood  estimation.  A linear,  time-invariant  analytical 
model  is  derived  in  detail.  Next,  the  cost  function  used  to 
evaluate  the  estimation  is  discussed.  Finally  the  Modified 
Ne wton-Raphson  minimization  algorithm  is  described. 

Chapter  3 discusses  data  acguistion  and  processing  as 
well  as  preliminary  inputs  to  the  estimation  program.  The 
aircraft  instrumentation  is  described  in  detail,  and  the 
steps  in  data  processing  are  outlined.  Specifying  aircraft 
dimensions  allows  nondimensionalization  of  the  estimated 
derivatives  for  more  accurate  comparisons.  Initial 
estimates  of  the  parameters  to  be  identified  are  important 
to  the  success  of  the  minimization  algorithm. 

Chapters  U and  5 detail  the  results  of 
lateral-directional  and  longitudinal  derivatives 
determination.  Initial  estimates  illustrate  problems  with 
the  data  and  the  analytical  model.  Steps  are  then  taken  to 
improve  the  quality  of  the  estimates.  Additional  runs  ire 
made  to  add  insight  into  problems  experienced. 

Chapter  6 summarizes  the  results  of  this  research  and 
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the  recommendations  for  further  research. 

Appendices  A and  B tabulate  the  results  of  runs 
described  in  Chapters  4 and  5 and  list  several  important 
averages.  Appendix  C lists  the  nomenclature  used  in  this 
thesis . 
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2.  THEORY  OF  PARAMETER  ESTIMATION 


Identifying  aircraft  parameters  from  test  flights 
begins  with  three  determinations: 

. Selection  of  the  analytical  model. 

. Identification  of  the  cost  function. 

. Determination  of  an  algorithm  to  minimize  the  cost 
function. 

Choosing  the  appropriate  analytical  model,  and  thus,  the 
parameters  to  identify  in  the  eguations  of  notion,  »ay 
involve  reducing  the  order  of  the  model  to  a nanagable  level 
and  making  simplifying  assumptions  to  reduce  the  complexity 
of  the  model.  A quantifiable  index  of  performance,  or  cost 
function,  permits  optimization  of  the  identification  process 
and  evaluation  of  the  success  of  the  process.  Maximum 
likelihood  estimation  uses  the  log  likelihood  function  of 
the  parameters  to  be  identified  as  the  basis  for  its  cost 
function.  A numerical  algorithm  which  reduces  the  cost 
function,  allows  digital  computation  of  the  parameter 
estimates  with  limited  and/or  no  need  for  operator 
judgement.  The  Modified  Newton-Raphson  algorithm  performs 
these  functions.  Figure  1 shows  the  relationship  of  these 
factors  to  the  process  of  parameter  estimation.  Each  of 
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these  factors  is  discussed  below. 

2 . 1 Equations  of  Motion 

Beginning  with  12  state  variables,  a nonlinear, 
time- varying,  6-degree-of-f reedom  model  describes  the  motion 
of  an  aircraft  in  an  inertial  reference  frame.  The  earth  is 
assumed  to  be  fixed  in  space,  and  earth  curvature  is 
neglected.  Variations  in  atmospheric  properties  and  gravity 
with  altitude  and  position  perturbations  of  interest  also 
are  assumed  negligible.  Aircraft  mass,  and  mass 
distribution  are  assumed  constant,  and  the  aircraft  is 
assumed  to  be  a rigid  body,  leading  to  the  nonlinear 
differential  equation: 

I (t)  = f (C,  x,u,w)  (1  ) 

where  f is  a vector  of  specific  forces,  moments  and 
kinematic  effects,  c is  a vector  of  aircraft  parameters  to 
be  identified,  x is  the  state  vector,  u is  the  control 
vector,  and  w is  the  disturbance  vector. 

Assuming,  in  addition,  that  the  aircraft  is  in  steady 
equilibrium  at  its  nominal  flight  condition  implies  that  the 
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The  6 dynamic 


equations  of  motion  are  time  invariant, 
equations  are  converted  from  the  inertial-axis  system  to  the 
body-axis  system  using  inertial-to- body-axis  transf or aation 
matrices  from  Ref  5.  Stochastic  inputs  , w , to  the  model 
are  assumed  negligible. 

The  equations  are  linearized  by  Taylor  series 
expansion,  taking  only  the  perturbation  solutions.  The  h 
state  variables  representing  "pure”  integrals  of  other 
varianles,  ( x y z ^ ) are  dropped  next,  since  they  add 
little  information  to  the  model.  Products  and  squares  of 
perturbation  quantities  are  assumed  negligible.  Sines  and 
cosines  of  perturbation  jangles  are  assumed  to  be 
approximately  the  angles  and  one,  respecti vely . The  linear 
differential  equation  model  derived  from  Eq  1 then  can  be 
expressed  as 


" al 

a2  ' 

x + 

H 

u 

(2) 

a3 

a 4 

L b2  J 

where  al  and  bl  are  longitudinal  mode  derivative  matrices, 
a4  and  b2  are  lateral-directional  mode  derivative  matrices,, 
and,  a2  and  al  are  mode  coupling  derivative  matrices. 

T 

x = Cuwgevpr*3  (3) 
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u = I 6e  ea  er  6sf  1 | 


Inertial  cross-products  between  the  longitudinal  and 

lateral  axes  and  the  lateral  and  normal  axes  are  assumed  to 

be  negligible  (I  = I = 0 ) . Straight  and  level, 

XY  YZ 

trimmed  flight  further  simplifies  the  equations,  since 


f?  = p = q = r=  * = 0 


(5) 


The  translational  states  are  converted  to  more  easily 


measurable  angles 

by 

the 

rela  tionshi ps: 

v / 

<u2. 

v2  + 

w2) 

1/2 

sin 

a 

u / 

(u2  ♦ 

V2. 

.2i 

1/2 

cos 

a 

(6) 

V / 

c2. 

v2  + 

»2) 

1/2 

sin 

0 

Figure  2 defines  the  body-axis-neasurable  aircraft  states 
and  control  deflections  in  their  positive  sense,  used  in  the 
equations  of  motion.  The  following  body-axis,  linear, 
t ime- invariant,  perturbation  differential  equations  are 
uncoupled  into  a longitudinal  motion  set  and  a lateral- 


2 


5 


directional  set 


u = [ «e  1 ] 


For  the  lateral-directional  case: 


Y *sina  Y -cosa  g/v  cos* 
P r 


cos*tane  0 


x = C 0 p r *3 
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u 


(15) 


u = [ 6a  6r  6sf 


T 

D 


where  A is  the  fundamental  matrix,  and  B is  the  control 
effect  matrix.  Dividing  the  equations  into  two  sets  still 
allows  determination  of  parameters  coupling  the  two  sets  as 
suggested  in  Pef  6.  The  division  reduces  computation  time, 
model  complexity  and  order,  simplifies  input  design,  and  is 
common  practice. 

Inclusion  of  certain  unsteady  aerodynamic  parameters, 

such  as  n , is  accomplished  using  the  acceleration 
a 

transformation  matrix,  K,  in  the  dynamic  equations: 

Px  = Ax  + Bu  (16) 


The  P.  matrix  is  the  identity  matrix  minus  the  matrix  of 
unsteady  aerodynamic  parameters. 

The  aircraft  analytical  model  is  made  up  of  three  sets 
of  linear  equations.  The  dynamic  equations  or  equations  of 
motion  (17)  were  just  derived.  The  output  equations  (18) 
define  the  observation  relationship  of  states  and  state 
rates.  The  observation  equations  (19)  incorporate  outputs 
and  measurement  noise.  Modeling  errors  and  instrument 
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biases  are  assumed  negligible 


(17) 

(18) 
(19) 

[ 

The  G and  H matrices  incorporate  the  A and  B natrices  and 
the  identity  matrix  to  produce  states  and  state  rates  in  the 
output  vector.  n is  the  measurement  noise  vector. 
Measurement  noise  is  assumed  to  be  a stationary,  Gaussian, 
white  noise  process.  Control  measurement  noise  is  assuted 
to  be  negligible. 

In  summary  , the  following  assumptions  are  observed  in 
the  aircraft  analytical  model: 

. Turbulence,  modeling  errors,  and  instrument  biases 
are  negligible. 

. Airframe  is  a rigid  body;;  i.e.,  there  are  no 
significant  structural  or  aeroelastic  nodes. 

. Earth  is  fixed  in  space. 

. Mass  and  mass  distribution  of  the  aircraft  are 
constant . 

. Products  and  squares  of  perturbation  quantities  are 
assumed  negligible. 
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. Variations  in  atmospheric  properties  with  altitude 

perturbations  of  interest  are  negligible. 

. Trim  values  are  p=p  = r = q = 4 = 0 

. The  inertia  cross-products  I and  I are 

XY  YZ 

negligible. 

. Measurement  noise  is  a stationary,  Gaussian,  white 
noise  process. 

. No  stochastic  inputs  are  allowed.  Control 
measurement  noise  is  negligible. 

2.2  Cost  Function 

The  objective  of  maximum  likelihood  estimation  is  to 

maximize  the  probability  density  function,  p(c|Z  ),  of  the 

N 

unknown  parameters,  c,  (stability  and  control  derivatives, 

state  initial  conditions,  aerodynamic  biases)  conditioned  on 

the  state  and  control  measurements,  Z . In  other  words,  the 

N 

estimator  must  maximize  the  likelihood  that  an  estimate  of 
the  "unknown  stability  and  control  derivatives  in  the 
equations  of  motion  will  produce  the  measured  state  time 


histories . 

Bayes's 

rule  expresses 

this 

conditional 

probability 

density 

function  in  terms 

of 

more  readily 

computed  quantities: 
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P (£|Z  ) = p(Z  |c)  p(c)  / p (Z  ) 
N N N 


(20) 


where  p (Z  |c)  is  the  probability  density  of  the  measurements 
N 

conditioned  on  the  unknown  parameters,  p (c)  is  the  prior 

probability  density  of  the  parameters,  and  p(Z  ) is  the 

N 

probability  of  the  measurements  and  equals  1. 

A sequential  computation  of  p(c|Z  ) can  be  derived  from 

N 

Bayes's  rule.  p (c)  is  approximated  by  p(c|Z  ),  the 

N - 1 

probability  density  of  the  parameters  conditioned  on  N- 1 
sample  time  measurements,  and  each  prior  conditional 
probanility  density  function  is  defined  in  a like  manner, 
back  to  the  starting  point  of  the  data.  Then, 


N 

P (c|Z  ) =TT  P(Z  |C)  p(c  1 0)  (21) 

N i = 1 N 


where  p(c|0)  -is  the  probability  density  of  the  parameters 
prior  to  any  measurements. 


If  measurement  errors,  process  errors,  ana  state 
initial  conditions  have  stationary,  Gaussian  probability 
distributions,  then  the  probability  distribution  of  the 


parameters  conditioned  on 

the 

measurements 

will 

also 

be 

Gaussian.  In  addition. 

the 

probability 

density 

of 

the 

measurements  conditioned 

on  the 

parameters  is 

the 

same 

as 
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the  probability  density  of  the  measurements  conditioned  on 
the  residual  calculated  using  the  parameters  (Bef  7)  : 

I T 

P(z  II  ) = M/2  -1  1/2  exp  r -1  (z  -I  ) D1  (z  -i  ) 1 

i i (2«)  | D 1 | L 2 i i i i J 

(22) 

where  D 1 is  the  measurement  error  covariance  matrix,  and  z 
is  the  observation  vector  at  sampling  time  i.  Taking 
advantage  of  the  properties  of  logarithms,  the  log 
likelihood  function  is: 

N T 

In  p (c  1 2 ) * -1  1 (z  ) Dl(z  ~J  ) 

N 2 i=1  i i i i 

N N/2  -1  1/2 

•J  2 In  (2*J  1 D 1 | 

2 i=  1 

-In  p (c ) 0)  (23) 

Por  a constant  error  covariance  natrix  and  a constant 
estimate  of  the  probability  density  of  the  parameters  prior 
to  any  measurements,  the  quadratic  cost  function,  J , which 
neglects  the  latter  two  terms  is  minimized  when  the  log 
likelihood  function  is  naxinized: 
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N T 

J = 1 $ (Z.  "I  ) D1  (Z  "I  ) 

N- 1 i-1  i i i i 


(2U) 


where  N is  the  number  of  sampling  instants, i is  the  sampling 
index,  and  D1  is  the  measurement  error  weighting  matrix 
(ideally  it  is  the  inverse  of  the  error  covariance  matrix) . 

This  mean-sguared-err or  criterion  serves  many  useful 
ends.  The  cost  function,  J,  compares  predicted  aircraft 
motion  with  measured  flight  time  histories.  It  reflects  the 
confidence  in  the  measurement  of  those  time  histories.  It 
weighs  large  errors  more  heavily  than  small  errors, and  it  is 
simple  to  use. 

There  are  several  ways  of  determining  the  D1  weighting 
matrix.  One  approach  is  to  base  the  weighting  on  published 
statistical  properties  of  the  instrument  noise.  Reference  7 
employs  a Kalman  Filter  state  estimator  to  estimate  the 
measurement  error  covariance  matrix,  then  inverts  it  to  form 
the  equivalent  of  D1.  Reference  3 and  8 suggest  using  the 
smallest  mea n-sguared-error  obtainable  from  the [test  flight 
data  as  a measure  of  the  noise  in  a measured  variable. 

In  this  study,  D1  matrix  elements  were  determined 
following  the  recommendations  in  Ref  6.  The  diagonal 
elements  of  D1  first  were  estimated  from  published 
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instrument  accuracy  and  resolution  information.  Next,  the 
cost  function  was  minimized  with  the  initial  Dl  matrix. 
Elements  of  the  Dl  diagonal  were  replaced  by  one  over  the 
mean-squared-error  in  a state  variable  if  the  weighted  mean 
squared-error  was  not  approximately  egual  to  one: 


d = 1 / a 
ii  ii 


i — 1 r • • • r 1 


(25) 


wnere  1 is  the  number  of  elements  in  the  observation  vector, 
z.  This  process  was  repeated  with  the  updated  Dl  matrix 
until  all  diagonal  Dl  elements  resulted  in  weighted  mean 
sq uared  — errors  near  one.  All  off-diagonal  Dl  elements  were 
set  to  zero  since  correlation  in  instrument  noise  was 
assumed  to  be  small. 

The  cost  function  was  modified  to  consider  ^ priori 
estimates  of  the  parameters  in  the  cost  function: 


N T 

J ~ _L  5 (2  ~1  ) D 1(z  "I  ) 

K- 1 i-1  i i i i 


(C-c  ) K D2  (c-c  ) 


(26) 


where  c is  a vector  of  all  parameters  to  be  identified,  c 


is 


the  vector  with  X priori  values  of  the  parameters  to  be 
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identified,  D2  is  a symmetric  weighting  matrix  reflecting 
the  relative  confidence  in  a priori  estimates,  and  K is  a 
constant  reflecting  the  importance  of  a departure  from  ^ 
priori  values  over  estimation  from  flight  data. 

A priori  weighting  permits  the  use  of  additional 
information  in  the  estimation  process.  Reference  9 suggests 
that  the  use  of  a priori  weighting  prevents  estimates  from 
changinq  from  a priori  values  unless  there  is  sufficient 
information  in  the  flight  data  to  justify  a change. 
Reference  7«  discusses  advantages  in  numerical  computations 
through  the  use  of  cl  priori  weighting.  The  cost  function, 
J,  as  defined  in  Eg  2b  forms  the  basis  of  the  Modified 
Maximum  Likelihood  Estimation  algorithm  which  is  used  here. 

2 . 3 Modified  Newton-Raphson  Minimization  Algorithm 

Reference  8 describes  several  parameter  identification 
methods  and  their  associated  cost  functions.  The  advantages 
of  using  digital  computation  to  reduce  the  cost  function  are 
also  ildescribed.  Digital  computation  more  readily  analyses 
large  amounts  of  data  with  little  or  no  operator  judgement 
reguired.  The  Modified  Newton-Raphson  minimization 
algorithm  approximates  the  first  and  second  gradients  of  the 
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cost  function  with  respect  to  the  unknown  parameters,  and 
from  these,  it  estimates  the  increments  in  the  parameters 
which  drive  the  cost  function  to  a minimum. 

The  Newton-Ra phson  iteration  technique  finds  a zero  of 
a nonlinear  function  of  several  parameters  (Ref  U)  . Setting 
the  gradient  of  the  cost  function  to  zero  defines  the 
minimum  cost  set  of  parameters.  Therefore, the 
Newton-Raphson  technique  can  be  used  to  derive  a 
minimization  algorithm. 

The  gradient  of  t^e  cost  function  can  be  expanded  in 
Taylor's  series  about  the  kth  iteration  value  of  c: 


2 

(V  J)  = (V  J)  + (V  J)  (c  -c  ) (27) 

c k ♦ 1 ck  ckk  + lk 

♦ Higher  Order  Terms 


where  (7  J)  is  the  first  gradient  of  the  cost  function 
c k 

with  respect  to  the  parameters  at  the  kth  .-.iteration,  and 

2 

(7  J)  " is  the  second  gradient  of  the  cost  function, 
c k 

Neglecting  higher  order  terms,  the  increment  of  the 
parameters  that  makes  (V  J)  approximately  equal  to  zero 
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(C  -c  ) 
k+ 1 k 


= -\  (V^J)  1 \v  J) 

L c k J c k ■ 


This  is  the  familiar  form  of  the  Newton-Kaphson  algorithm. 
The  first  and  second  gradients  of  the  cost  function  are 
defined  in  terms  of  the  gradients  of  the  residual: 


N T 

(7  J)  = 2 ^ (1  "I  ) D1  V (z  -j£  ) 

c k N-1  i = 1 i i k c i i k 

2 N T 

(V  J)  = _2_  ^ V (z  -i  ) D1  7 (z  -j  ) 
c k N-1  i=1  ciik  ciik 

N T 2 

+ _2_  ^ (z  'j;  ) D1  7 (z  -j  ) 

N- 1 i = 1 i ik  c i ik 


with 


(z  - I)  = (2  - 1)  * 7 (2  - 2)  (C  - C ) 

k k-1  c k k+1 


where  V (z  -y  ) is  the  first  gradient  of  the  residual  for 
ciik 

the  kth  iteration  and  the  ith  sampling  instant,  and 

2 

V (z  -jr  ) is  the  second  gradient  of  the  residual  tor  the 
ciik 

kth  iteration  and  the  ith  sampling  instant. 

Calculation  of  the  first  gradient  of  the  residual  is 
straightforward  (Ref  8)  and  is  analogous  to  the  calculation 
of  the  "sensitivity  equations"  in  Ref  10: 
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7 (I  - I)  = - (V  G)  x - G (V  x)  - (7  H ) u -H  (V  u) 
c c c c c 

(32) 

Balakrishnan  shows  in  Ref  11  that  the  second  gradient 
of  the  residual  goes  to  zero  as  the  parameter  estimates 
approach  their  true  values.  Neglecting  this  term  reduces 
the  amount  of  computation  required,  and  still  results  in 
asymptotically  unbiased  estimates.  The  algorithm  is  now 
referred  to  as  the  Hodified  Newton-Raphson  algorithm. 

Reference  8 discusses  the  use  of  an  approximation  to 
the  Cramfer-Rao  bound  in  conjunction  with  rthe  Modified 
Newton-Raphson  minimization  algorithm  used  here  to  indicate 
confidence  levels  for  the  parameter  estimates.  If  the 
actual  parameter  covariance  matrix  were  known,  it  would 
indicate  which  parameters  had  been  most  reliably  estimated. 
Reference  11  shows  that  this  covariance  matrix  is  bounded 
from  below  by  the  matrix  Cram^r-Rao  bound.  By  assuming  that 
the  estimates  obtained  using  the  Modified  Newton-Raphson 
algorithm  are  asymptotically  unbiased  and  that  the 
measurement  noise,  n,  is  a stationary,  Gaussian,  white  noise 
process,  the  matrix  Cram^r-fiao  bound  is  the  inverse  to  the 
second  gradient  of  the  cost  function  matrix  with  respect  to 
the  parameters.  Similar  assumptions  are  used  in  Ref  10  to 
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obtain  an  approximation  to  the  parameter  covariance  matrix. 
The  Cram^r-Sao  lower  bounds  to  the  standard  deviations  in 
the  parameter  estimates  are  used  in  this  thesis  to  indicate 
the  confidence  in  estimates  from  different  sets  of  time 

histories. 

Reference  4 is  a listing  of  a program  written  in 

FORTRAN  IV  using  the  equations  of  motion,  cost  function,  and 

minimization  algorithm  here  described.  Authored  by  Maine 
and  Iliff,  it  has  been  used  to  analyze  thousands  of 

maneuvers  on  several  aircraft. 

This  program  generates  state  time  history  estimates 
using  the  state  transition  matrix  defined  in  Ref  12: 

x(t  ) = * X(t  ) ♦ f +(t  ,T)  B U (t)  dT  (33) 

2 1^1 

where  * is  the  state  transition  matrix  and  u is  the  measured 
control  vector. 

Normal  and  lateral  acceleration  time  history  estimates 
are  made  using  Eq  18.  The  specific  equations  are  restated 
here: 

a = -VZor-VZ  Be  - V Z ♦ cos#,  cose  (34) 

n 9 a g Be  g o 
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a = V Y £ ♦ V Y ea+VY  er  4 V Y 6sf  4 V Y (35) 

y 9 /J  q ea  g 6r  g esf  g o 


Control  leasureiient  noise  is  assumed  to  be  negligible. 
The  validity  of  this  assumption  is  evident  in  the  guality  of 
the  parameter  estimates. 

This  program  was  used  with  only  minor  modification  on 
VRA  flight  data  collected  from  two  test  flights.  The  next 
chapter  describes  the  collection  and  processing  of  that  data 
and  the  preliminary  inputs  to  the  computer  program. 
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3.  DATA  ACQOISITION  AND  PRELIMINARY  INPUTS  FOK  PARAMETER 


ESTIMATION'. 

Princeton’s  VBA  has  a full  complement  of  instruments 
and  necessary  telemetry  to  measure  all  states  and  some  state 
rates  of  interest.  Upon  calibrating  this  instrument 
package,  flight  data  were  collected  and  recorded  for 
processing.  A nalog- to-digita 1 conversion,  application  of 
calibrations,  and  correction  of  time  shifts  prepared  the  raw 
data  for  use  in  the  estimation  program.  In  addition,  true 
airspeed  and  dynamic  pressure  were  calculated,  and  time 
histories  were  selected. 

Specifying  aircraft  dimensions  allowed  nondimensional- 
ization  of  derivatives  for  comparison  with  published  values. 
The  measurement  error  weighting  matrix  was  approximated 
using  several  time  histories.  Preliminary  estimates  of  the 
unknown  parameters  were  taken  from  unpublishej  analog 
matchinq  results  and  published  estimates. 

3 . 1 Instrumentation 

In  the  aircraft,  instrument  readings  and  control 
surface  positions  are  sampled  at  twenty  and  forty  samples 
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The  VRA  has  two  an gle-of -attack  vanes,  ODe  aounted  on  a 
boom  extendinq  from  each  wing,  at  the  quarter  chord.  This 
position  does  not  require  an  angular  rate  correction,  but 
previous  test  flights  indicate  that  an  aerodynamic  scale 
factor  of  1.6  must  be  used  to  correct  for  upwash  effects. 
The  signals  of  the  two  vanes  are  averaged  prior  to 
transmission  via  telemetry. 

Perturbations  in  the  wing  vortex  flow  with  changing 
angle-of -attack  and  aileron  deflection  caused  a dynamic 
error  in  the  single  sideslip  i angle  v ane  mounted  on  the 
right  wing  boom.  In  Figure  3,  the  positive  sideslip  anqle 
perturbation  (H)  of  approximately  1-deg  at  1-sec  is  nearly 
in  ^.phase  with  the  aileron  pulse  input.  The  sideslip  angle 
estimate  (E)  based  on  present  study  stability  and  control 
derivatives  does  not  predict  the  measured  perturbation. 
Examination  of  the  lateral  acceleration  measurement  shows 
some  acceleration  due  to  aileron  deflection  but  it  is 
questionable  whether  the  magnitudes  of  the  perturbations  are 
in  agreement.  Since  the  aileron  pulse  time  histories  came 
from  the  first  test  flight,  side-force  derivative  estimates 
using  those  time  histories  are  suspect.  After  the  first 
flight  an  additional  vane  was  added  to  the  left  wing  boom, 
and  both  signals  then  were  averaged  prior  to  telemetry 
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FIGURE  3.  COMPARISON  OF  SIDESLIP  ANGLE  AND  LATERAL 
ACCELERATION  MEASUREMENTS  (M)  FOR  A FLIGHT  1 AILERON  INPUT 
TIME  HISTORY  (E-  STATE  ESTIMATES  BASED  ON  PRESENT  STUDY 
STABILITY  AND  CONTROL  DERIVATIVES).  (ABSCISSA  IS  TIME  IN 
SEC) 
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transmission.  This  correction  must  be  evaluated.  An 
anqular  rate  correction  for  the  vanes'  1.3  ft  displacement 
above  c.g.  is  employed. 

The  attitude  gyros,  pitch  rate  gyros,  and  normal  and 
lateral  accelerometers  are  clustered  about  the  main  wing 
spar,  on  or  near  the  c.g.  This  location  reduces  the 
sensitivity  of  the  pitch  rate  gyro  to  structural  vibration. 
On  Flight  1,  a single  accelerometer  was  available  and  it 
measured  lateral  accelerations.  A second  accelerometer  was 
installed  prior  to  Flight  2.  Figure  4 ;shows  a typical 
lateral  acceleration  time  history  from  Flight  2.  The  .25  g 
noise  band  in  the  measurement  (H)  is  not  apparent  for  the 
Flight  1 lateral  acceleration  time  history  in  Figure  3. 
Since  the  lateral  acceleration  measurement  was  not  useful, 
tne  guality  in  the  rDutch  roll  related  stability  and  -control 
derivatives  was  degraded.  Further  investigation  into  the 
causes  of  this  noise  are  required. 

The  roll  and  yaw  rate  r gyros  are  mounted  in  the 
equipment  bay  behind  the  c.g.  The  poor  lateral  acceleration 
signal  and  the  failure  of  the  roll  rate  gyro  degraded  the 
estimation  of  side-force  and  rolling  moment  derivatives  from 
Flight  2 data.  Estimation  of  lateral  derivatives  still  was 
attempted  since  sideslip  angle  and  bank  angle  time  histories 
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FIGURE  4.  NOISY  LATERAL  ACCELERATION  MEASUREMENT  (H)  FROM 
FLIGHT  2 SIDE-FORCE  INPUT  TIME  HISTORY  (E-  STATE  ESTIMATES 
BASED  ON  PRESENT  STUDY  STABILITY  AND  CONTROL  DERIVATIVES)  . 
(ABSCISSA  IS  TIME  IN  SEC) 


were  available. 

The  aircraft  pitot  boon  is  mounted  below  and  forward  of 
the  right  wing*  The  two  static  ports  are  on  either  side  of 
the  fuselaqe.  Large  pitch  rates  accompanying  elevator 
doublets  produced  a noticeable  instrument  response  in  the 
velocity  measurement  as  seen  in  Figure  5.  The  large  , 
velocity  perturbation  starting  at  2-sec  corresponds  to  a • 2 
g deceleration  that  was  not  confirmed  by  the  pilot*  This 
dynamic  velocity  error  is  produced  by  the  changing  pressure 
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measured  by  the  static  port  with  rapid  changes  in  flow  over 
the  wing.  Although  this  problem  is  not  a serious  concern  to 
production  models  of  the  Navion,  the  dynamic  error  prevented 
use  of  the  velocity  time  history  in  estimating  longitudinal 
derivatives  from  elevator  doublets. 

3.2  Calibration 


Aircraft  instruments  and  telemetry  signals  were 
calibrated  prior  to  each  flight.  Deflecting  the  instruments 
and  control  surfaces  to  measured  positions  or  rates,  tne 
telemetry  siqnals  were  recorded  in  units  of  percent 
telemetry.  These  units  correspond  to  the  percentage  of  the 
maximum  pulse  duration  transmitted  via  telemetry.  Scale  .and 
bias  factors  were  calculated  for  each  instrument  and  control 
surface  using  linear  regression.  Table  2 lists  signal  scale 
and  bias  factors  for  Flights  1 and  2.  The  angle-of-attack 
factors  reflect  corrections  for  upwash  effects. 

Overall,  there  is  little  change  in  scale  factors  and 
biases  between  the  two  flights,  a surprising  result 
considering  that  the  instruments  were  replaced  during  the 
3-month  period  between  flights,  that  there  was  normal  wear 
and  tear  of  the  data  acquisition  system  with  usage,  and  that 
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a number  of  people  were  involved  in  making  the  calibrations. 

The  bias  differences  for  the  angle-of -attack  and 
sideslip  vanes  were  expected  due  to  differences  ifc  alignment 
procedure.  The  vanes  were  faired  with  the  wing  for  the  zero 
reference  values  prior  to  Flight  1,  and  they  aligned  to 
produce  17  "units"  on  the  angle-of-a ttack  meter  and  zero 
"units"  on  the  sideslip  meter  as  zero  reference  values  prior 
to  Flight  2. 

The  scale  factor  and  bias  differences  for  the  bank 
angle  attitude  gyro  are  of  concern.  Only  4 data  points  were 
used  to  calculate  the  scale  factor  and  bias  for  the  first 
test  flight,  and  it  is  possible  that  2 data  points  were 
beyond  the  linear  ranqe  of  the  ground  station.  Since  the 
response  of  the  attitude  gyro  for  Flight  2 was  not  similar 
to  the  data  for  Flight  1 , it  was  not  felt  that  using  the 
scale  factor  and  bias  from  Flight  2 was  superior  to  using 
questionable  values  from  Flight  1.  Comparison  of  roll  rate 
and  roll  angle  time  histories  for  Fliqht  1 did  not  point  to 
any  gross  discrepancies. 

The  lateral  accelerometer  values  represent  two 
different  accelerometers,  telemetered  on  two  different 
channels.  In  light  of  the  tact  that  the  lateral 
acceleration  time  histories  for  Fliqht  2 are  very  noisy. 
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TABLE  2.  CALIBRATION  SCALE  FACTORS  AND  BIASES  FOR 
FLIGHTS  1 AND  2. 


Scale  Factor  Bias 

Percent 

Fit.  1 Fit.  2 Change  Fit.  1 

Fit.  2 

Percent 

Change 

a 

-.1362 

-.1315 

2 

8. 30 

5.88 

17 

q 

-.4725 

-.4892 

2 

20.77 

22.  19 

3 

V 

-.7031 

- 

129.79 

- 

e 

-.3  977 

-.3917 

1 

19.30 

20.43 

3 

a 

.0153 

- 

-.7970 

— 

n 

P 

-.4385 

-.4301 

1 

20.85 

17.49 

9 

P 

-.5318 

-.5473 

1 

24.  34 

24.34 

2 

r 

.4290 

.4399 

1 

-19.38 

-20.45 

3 

<*> 

.7577' 

.8934 

8 

-37.36 

-40.87 

4 

a 

.0170 

.0176 

2 

-.7273 

-.8196 

6 

Y 

66 

-1.0220 

-.9150 

5 

38.11 

31.99 

9 

ea 

- . 5295 

-.4862 

4 

23.29 

21.24 

5 

6r 

. 3750 

.4396 

8 

- 19.  19 

-25.76 

15 

6 S f 

-.8101 

-.7991 

1 

43.50 

44.42 

1 

pointing  to  some  installation  error  or  hardware  problea,  it 
is  surprising  that  these  tens  are  as  close  as  they  are. 

The  relatively  large  difference  in  control  surface 


I 


scale  factors  and  biases  between  flights  was  not  expected. 
No  new  potentiometers  and  no  extensive  rewiring  occurred 


between  flights.  Since  the  state  and  state  rate  factors 
changed  little,  it  is  unlikely  that  the  telemetry  package  or 
the  qround  station  produced  the  differences.  The  side-force 
surrace  calibrations  used  protractor  markings  permanently 
etched  on  the  wing  which  would  tend  to  reduce  human  error  in 
measurement.  The  elevator  and  aileron  deflection  angles 
were  measured  using  hand-held  bubble  protractors,  and  the 
rudder  measurements  for  Flight  1 were  made  using  a damaged 
protractor.  Human  error  in  measuring  surface  deflections 
might  be  reduced  by  using  protractors  specially  constructed 
for  each  surface  and  bolted  in  place  during  calibrations. 
Such  a protractor  has  already  been  constructed  and  is  now 
used  for  the  rudder  deflection  measurements.  The  ease  of 
using  this  protractor  suggests  that  human  error  is  being 
reduced . 

3 . 3 Data  Collection 

The  data  used  in  this  study  were  obtained  from  two  test 
flights  of  the  VFA  at  105  KIAS.  Outside  air  temperature  and 
pressure  altitude  were  noted  for  later  use  in  calculating 
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true  airspeed  and  dynamic  pressure.  The  applied  control 
inputs  included  elevator  doublets  and  pulses,  aileron 
pulses,  rudder  doublets,  and  side-force  surface  pulses  and 
doublets. 

During  Flight  2,  aircraft  pitch  angle  and 
angle-of-atta ck  were  measured  using  a bubble  protractor  and 
angle-of-attack  meter.  Both  indicated  a zero  pitch  angle 
and  angle-of-attack  (within  1/2  deg)  with  respect  to  the 
fuselage  reference  line  for  straight-and-level  flight  at  105 
KIAS. 

The  aircraft  nominal  altitude  was  selected  during 
flight  to  minimize  turbulence.  105  KIAS  was  the  flight 
airspeed  used  for  two  other  research  programs  conducted 
concurrently  with  this  study.  It  guarantees  low  trim 
angle-of-attack  to  remain  within  the  linear  range  of  the 
lift  curve  slope  during  control  inputs.  In  addition,  the 
control  surfaces  are  more  effective  than  at  slower  speeds, 
thus  producing  state  perturbations  with  good  signal-to-noise 
ratios  (Ref  13).  The  guality  of  the  control  inputs  reflects 
the  experience  of  the  test  pilot  with  analog  matching 
parameter  estimation. 

The  VBA  is  now  eguipped  with  a Microprocessor  Digital 
Flight  Control  System  that  has  demonstrated  its  ability  to 
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produce  precise  square  elevator  pulses.  For  future 
parameter  estimation  studies,  thi.s  system  could  be  used  to 
input  precise  control  deflections  better  suited  to  excite 
several  aircraft  modes  of  motion  simultaneously  or 
sequentially  (Ref  14).  In  addition  optimal  input  desigD 
studies  could  be  demonstrated  using  the  VRA  (Ref  3,16). 

% 

3 . 4 Analoq-to-Diqital  Conversion 

Flight  test  data  were  recorded  during  flight  on  a 
Honeywell  76HC  1-inch  tape  recorder.  Due  to  equipment 
compatibility  requirements  the  data  were  re-recorded  on 
1/4-iDCh  tape.  The  tape  was  transferred  to  another  facility 
to  be  digitized.  The  discrete  analog  signals  on  the  tape 
were  converted  to  "boxcar"  analog  siqnals  using  integrator 
circuits  in  the  ground  station  (Ref  16).  The  signals  were 
sampled  simultaneously  at  4C  SPS  per  signal  and  digitized  by 
a SSL  600  Analog-to-Digita 1 Converter.  The  digitized 
samples  were  recorded  in  6-bit  BCD  on  7 track,  160C  bit- 
per-inch  magnetic  tape  usinq  an  AHPEX  Digital  Recorder. 

The  resulting  tape  cannot  be  processed  directly  by  a 
FORTRAN  program  because  PORTRAN  does  not  recognize  the 
manner  in  which  signs  of  samples  are  stored  on  the  tape. 
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The  "Seven  to  Nine  Track  Tape  Conversion  Program",  written 
by  Richard  B.  Gilbert,  Technical  Staff  Member  in  the 
Department  of  Mechanical  and  Aerospace  Engineering,  reads  the 
tape,  scales  the  data  to  percent,  telemetry  and  records  it  in 
an  acceptable  FORTRAN  format. 

This  analog-to-digital  conversion  procedure  is  time 
consuming  and  it  adds  several  sources  of  errors  in  the  raw 
data.  Each  step  in  the  process  was  delayed  by  failures  in 
aging  equipment.  Plans  are  being  made  to  replace  this 
equipment.  The  replacement  system  should  digitize  the 
discrete  analog  samples  directly  from  the  original  data 
tape.  This  would  prevent  the  addition  of  noise  to  the  data 
by  the  additional. tape  recorders  and  the. ground  station.  It 
also  would  eliminate  the  ground  station  as  a source  of  bias 
errors. 

3.5  Data  Scaling 

After  producing  a magnetic  tape  containing  digitized 
fliqhb'  time  histories  in  an  acceptable  FORTRAN  format  the 
data  must  be  converted  into  engineering  units,  interpolated 
to  reduce  the  errors  produced  by  time  shifts,  and  corrected 
for  compressibility  and  density  effects.  This  final  data 
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processing  was  performed  using  a FORTRAN  program  entitled 


♦Sea  ling  * . written  as  part  of  this  '•  research.  The  program 
first  applies  the  calibration  scale  factors  and  biases  to 
convert  the  digitized  data  from  percent  telemetry  to 
engineering  units. 

Scali nq  calculates  true  airspeed  and  dynamic  pressure 
from  measured  indicated  airspeed  and  pilot  judged  nominal 
pressure  altitude  and  outside  air  temperature.  The 
equations  used  in  these  calculations  were  derived  from  the 
perfect  gas  law  and  definitions  of  the  standard  atmosphere. 
They  are  taken  directly  from  Ref  17: 


2 

q = 1 r V 

2 true 


(37) 


where  Ta 
calibrated 
and,  p is 
divided  by 


is  absolute 
airspeed,  H 
air  density 
Ta  and  309D. 


c 


air 

is 

and 


temperature (deg  K) , V is 

c 

calibrated  pressure  altitude, 
equals  atmospheric  pressure 
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Reference  18  discusses  the  effects  of  tine  delays 
between  sampling  instants  of  states  and  controls  on  the 
performance  of  a maximum  likelihood  estimator.  Steers  and 
Iliff  concluded  that  delays  of  control  measurements  caused 
the  greatest  degradation  in  the  accuracy  of  estimates. 
Since  no  measured  quantities  were  sampled  simultaneously 
prior  to  transmission  to  the  ground,  and  control 
measurements  lagged  state  measurements  by  as  much  as  30 
percent  of  the  sampling  period,  it  was  deemed  important  to 
interpolate  between  sampling  instants  to  produce  data  with 
negligible  time  shifts.  Seal in  g also  performs  this 
interpolation  on  a measurement,  s,  using  the  following 
eq  ua  tion: 

s = s ♦ k (s  -s  ) (38) 
i i i ♦ 1 i 

where  k = 1 - time  delay/sampling  period  (39) 

The  analoq-to-digital  converter  samples  all  quantities 
simultaneously  and  therefore  does  not  add  additional  delays. 
Since  sideslip  angle  is  the  first  sampled  signal, all  data 
were  interpolated  to  the  time  corresponding  to  the  sideslip 
angle  sampling  instant. 
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Finally,  Scaling  calculates  means  and  standard 
deviations  of  all  states,  state  rates,  and  controls  measured 
as  an  aid  in  determining  biases.  These  biases  reflect  fixed 
biases  used  by  the  ground  station  and  analog-to-digital 
converter  in  manipulating  the  data,  fixed  instrument  biases 
in  fliqht  (for  example,  the  trim  angle-of -attack  vane 
position  is  fixed  by  flow  over  the  wing  in  trimmed  flight) , 
and  instrument  drift  between  calibration  of  the  instruments 
and  recording  of  the  data.  Fixed  biases  in  states,  state 
rates,  and  controls  are  inputs  to  the  Modified  Maximum 
Likelihood  Estimation  program. 

3 . 6 Inputs  To  The  Estimation  Program 

Time  histories  to  be  used  in  the  estimation  "program 
were  selected  primarily  by  two  criteria.  First,  the  time 
history  had  to  be  free  of  dropouts  and  other  eguipment 
produced  ' anomalies  as  discussed  in  Ref  9.  Second,  control 
inputs  were  chosen  that  appeared  to  be  greater  than  5 times 
the  noise  band  (peak  to  peak)  of  the  control  measurement  to 
insure  a good  signal-to-noise  ratio  (ratio  of  largest  signal 
perturbation  to  noise  band)  . 

Table  3 lists  aircraft  dimensions  and  flight  conditions 
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TABLE  3.  AIRCRAFT  DIMENSIONS  AND  FLIGHT  CONDITIONS 


A 


Hinq: 

Area 

180.0 

FT 2 

Span 

33.  38 

FT 

C hor  d 

5.  67 

FT 

Leading  Edge  Sweep 

3.0 

DEG 

Dihedral 

7.5 

DEG 

Hoot  Incidence 

2.0 

DEG 

Tip  Incidence 

Airfoil: 

-1.0 

DEG 

Tip 

NACA  641  OR 

Foot 

NACA  4415P 

Aileron: 

2 

Area 

5.4 

FT 

Deflection 

± 20.0 

DEG 

Plap: 

2 

Area 

83.6 

FT 

Deflection 

1 30.0 

DEG 

Horizontal  Tail: 

o 

Area 

43.0 

FT* 

Incidence 

-3.0 

DEG 

A irf oil 

NACA  0012 

Elevator: 

2 

Area 

14.1  . 

FT 

Deflection 

UP  30, DOWN  20 

DEG 

Vertical  Tail: 

2 

Area 

18.  1 

TT 

A irfoil 

MODIFIED  NACA  0011.2 

Fin  Offset 

2.0'- 

DEG 

Rudder : 

? 

Area 

11.6 

FT 

Deflection 

LEFT  23, RIGHT  17 

DEG 

Side  Force  Surface: 

2 

Area 

16.0 

FT 

Airfoil 

NACA  0012 

Deflection 

± 30.0 

DEG 

Gross  Height 

2609. 0 

LB  5 

I 

1573. 7 SLUG-FT 

X 

2 

I 

2736.0  SLUG-FT 

¥ 

I 

3673. 8 SLUG-FT2 

Z 

TABLE  3.  (CONTINUED) 


Flight  1 

Flight  2 

Pressure  Altitude 

3500  FT 

650u  FT 

Nominal  Indicated  Airspeed 

105  KT 

105  KT 

Outside  Air  Temperature 

0 degC 

10  degC 

Flap  Setting 

10  DEG 

10  DEG 

tor  Flights  1 and  2.  Tnis  information  is  used  to 

nondi mensiona lize  stability  and  control  derivatives. 

Side-force  surfaces  capable  of  producing  1/2  q lateral 

acceleration  for  full  deflection  at -.105  KIAS  were  added  in 

recent  years.  In  conjunction  with  this  modification,  the 

rudder  area  was  increased.  An  important  goal  of  this 

research  was  to  identify  the  side-force  surface  control 

derivatives.  Moments  of  inertia  were  modified  to  reflect 

the  weight  and  position  of  these  surfaces. 

No  actual  measurements  were  made  to  determine  the 

aircraft  gross  weight  and  moments  of  inertia.  Such  an 

experiment  would  be  of  great  benefit  in  light  of  the 

structural  modifications  made  to  the  aircraft.  In  addition, 

no  estimate  of  the  cross  moment  of  inertia  about  the 

longitudinal  and  lateral  axes, I , is  available  from 

xz 

previous  research.  I appears  in  several  elements  of  the 

XZ 

fundamental  matrix.  A,  and  although  an  1 of  zero 

XZ 


too 


simplifies  the  equations,  this  assumption  may  be  too 
restrictive  and  add  to  the  variance  in  estimating  these 
deri vatives. 

Although  theoretically,  the  minimization  algorithm 
should  converge  on  a set  of  parameter  estimates  given  zero 
values  of  the  initial  estimates,  convergence  on  the  true 
values  of  the  parameters  is  dependent  on  initializing  the 
estimation  program  with  a reasonably  good  set  of  ^ priori 
parameter  values  (Ref  3,7,8). 

Unpublished  stability  and  control  derivatives  which 
were  estimated  using  analog  matching  i techniques  are 
available  for  the  VPA  in  its  current  configuration. 
Reference  1,  2,  and  19  (quoted  in  Ref  3)  also  provide 
estimates  for  the  Navion  aircraft  in  various  configurations. 
The  analog  matching  results  proved  to  be  adequate  as 
starting  estimates. 

The  measurement  error  weighting  matrix  was  established 
for  each  set  of  control  input  time  histories  using  the 
method  described  in  Chapter  2. 

3 . 7 Conclusions 

Data  acquisition  is  a crucial  aspect  in  a parameter 


3 -20 


r 


estimation  process.  The  quality  of  the  flight  data  directly 
effects  the  quality  of  the  parameter  estimates.  Dynamic 
errors  in  the  sideslip  angle  and  velocity  measurement 
instruments  suggest  that  these  instruments  should  be 
repositioned.  Low  signal- to-noise  ratios  for  accelerometers 
suqgest  that  some  form  of  filtering  prior  to  transmission  of 
these  measurements  is  desirable.  In  calibrating  instruments 
and  control  surfaces,  specially  constructed  protractors  and 
care  in  avoiding  the  linear  limits  of  the  telemetry  and 
ground  station  will  reduce  scale  factor  and  bias  errors. 
Use  of  existing  equipment  on  board  the  VR A can  improve  the 
quality  of  control  inputs  for  use  in  parameter  estimation. 
The  discrete  analog  sampled  data  transmitted  to  the  ground 
should  be  digitized  directly  to  avoid  the  addition  of  noise 
and  biases  from  processing  data  throuqh  several  steps.  Time 
delays  introduced  by  the  telemetry  system  can  be  corrected 
by  interpolation  as  part  of  the  data  processing. 

Time  histories  selected  for  use  in  estimating 
derivatives  should  have  good  s ignal-to-noise  ratios  and  be 
free  of  equipment  related  anomalies.  Accurate  measurements 
of  aircraft  mass  and  inertia  properties  should  be  made  prior 
to  (or  as  part  of)  further  research  in  parameter  estimation. 
Published  estimates  and  previous  analog  matching  results 
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make  good  initial  estimates  for  the  estiBation  program. 

Chapters  4 and  5 discuss  tne  lateral-directional  and 
longitudinal  stacility  and  control  derivatives  estimated 
using  the  Modified  Maximum  Likelihood  Estimation  program. 
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U . LATERAL -PI SECTIONAL  DERIVATIVE  PETERBI NATION 

The  VRA  lateral-directional  derivatives  were  estimated 

using  twenty  time  histories  representing  four  basic  control 

inputs:  aileron  pulses,  rudder  doublets,  side-force  surface 

pulses,  and  side-force  surface  doublets.  The  measured 

quantities  were  sideslip  angle,  roll  rate,  yaw  rate,  bank 

angle,  lateral  acceleration,  aileron  deflection,  rudder 

deflection,  and  side-force  surface  deflection.  The  results 

of  parameter  estimation  runs  using  these  time  histories  were 

used  to  generate  statistics  describinq  the  quality  of  the 

estimates.  Initially,  the  variances  of  all  estimates  were 

high  v so  a number  of  procedures  were  use  to  determine  the 

quality  of  the  data  and  to  reduce  the  variance.  Results 
» 

were  separated  by  type  of  input  and  initial  direction  of  the 
primary  control  surface.  Because  the  variances  in  certain 
derivatives  were  reduced  by  separating  runs,  the  effects  of 
poor  siqna 1- to- noise  ratios  and  nonlinear  aircraft  responses 
were  documented. 

Four  initial  runs  indicated  deficiencies  in  the  data. 
Next,  aileron  input  time  histories  were  used  to  generate  a 
set  of  starting  values  for  the  rudder  time  history  estimates 
and  to  estimate  aileron  control  derivatives.  Rudder  input 
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time  histories  were  used  holding  aileron  control  derivatives 
fixed  since  the  signal-to-noise  ratio  of  aileron  deflections 
in  these  time  histories  was  small.  Using  the  aileron  and 
rudder  control  derivatives  already  estimated  from  previous 
runs,  estimates  were  made  using  the  side  force  surface  input 
time  histories.  Here  both  aileron  and  rudder  control 
derivatives  were  held  fixed.  From  these  runs,  an  overall 
set  of  estimates  emerged,  and  the  set  was  compared  with 
results  from  other  studies. 

h . 1 Initial  Huns 

Lateral-directional  parameter  estimation  began  with 
data  from  Flight  1.  Control  inputs  used  included  aileron 
pulses,  . rudder  doublets,  and  side-force  doublets. 
Deflection  of  the  side-force  surfaces  produced  rudder  and 
aileron  deflections  as  well,  due  to  a side-force-surf ace- to- 
rudder  interconnect,  and  a side-force-surf ace-to-aileron 
interconnect.  Parameter  estimation  using  each  type  of  input 
pointed  to  certain  deficiencies  in  the  data.  Using  a 5-sec 
time  history  of  a u-deg  aileron  pulse,  the  minimization 
algorithm  diverged.  The  minimization  for  a 6-sec  time 
history  of  a +2.5-deq  rudder  doublet  also  diverged.  The 
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siqnal-to-noise  ratios  for  these  control  inputs  were 
approximately  R and  5 for  a noise  band  or  1/2  deg  for  both. 
Larger  amplitude  inputs  resulted  in  the  minimization 
algorithm  converging  upon  a set  of  estimates.  The  estimates 
from  tour  time  histories,  two  aileron  pulses  ( + 7-deg  and 
-10.7-deg,  respectively),  a side-force  surface  pulse  of 
-3.5-deg,  and  a series  of  side-force  surface  doublets 
(^2.5-deg  each),  are  presented  in  Table  4 along  with  the 
standard  deviation  as  a percentage  of  the  estimate. 

TABLE  4. LATERAL  DIRECTIONAL  DERIVATIVE  ESTIMATES  AND  PERCENT 
STANDARD  DEVIATIONS  FOR  SEPARATE  CONTROL  SURFACE  INPUTS 


I ni tia  1 

Aileron 

Rudder 

Side  Force  Surface 

Runs 

Pulses 

Do  ublets 

Fit.  1 

Fit.  2 

C -.7878 

-.7795 

-. 875 8 

- 1.  002 

- 1.  155 

Yr  33.3 

25.1 

10.7 

33.9 

30.9 

C -.0812 

.0230 

.1667 

Ysa  206.8 

96.9 

267.2 

C .2043 

- . 0790 

. 3634 

. 2306 

J6 r 260.7 

293.5 

8.6 

127.7 

C .0861 

-.3252 

- . 1888 

.4499 

. 5837 

sf  305.8 

49.0 

229.3 

16.0 

14.7 

3 


TABLE  4.  (CONTINUED)  .' 


Initia  1 

A ileron 

Rudder 

Side  Force  Surface 

Runs 

P ulses 

Doublets 

Fit . 1 

Fit. 2 

C, 

-.0823 

-.0892 

-.1106 

-.1232 

.0877 

x0 

19.0 

16.0 

19.7 

14.5 

42.2 

C, 

-. 5012 

\ 

- .5609 

-.6376 

-.6858 

-.5562 

lp 

20.8 

8.8 

22.3 

21.4 

33.0 

C, 

.0700 

. 1103 

. 1088 

.0201 

. 1424 

1T 

199.6 

61.3 

15.8 

468.0 

96.2 

C- 

-.0709 

-.1818 

-.1272  /. 

*6a 

87.7 

13.1 

17.6  . 

C. 

- . 0472 

-.0379 

.0246 

. G24  6 

16r 

43.  1 

60.5 

5.2 

4.  1 

c. 

.0117 

. 0457 

. 1152 

.0069 

. 0040 

A6Sf 

27C.5 

93.9 

36.2 

121.2 

441.2  i. 

C 

.0875 

. 0745 

.0751 

.0751 

.0900  J 

P 

24.6 

18.4 

20.7 

15.3 

25.0 

c 

. 1051 

. 0865 

-.0935 

-.1327 

-.1755 

np 

268.4 

23.0 

53.0 

9.2 

76.2 

C 

- . 1310 

-.1211 

-. 1556 

-.2168 

-.0951 

nr 

14.5 

21.0 

14.6 

51.1 

101.8 

C 

-.0109 

-.0196 

-.0281 

nea 

62.3 

30.  1 

60.6 

c 

- . 0019 

-.0111 

- .0974 

-.0710 

n 

er 

21  12 

362.9 

4.3 

1C 2.  3 

c 

-.0020 

.0150 

- . 0029 

. 0430 

.0430 

6Sf 

768. 

125.7 

980..  0 

11.4 

17.1 

The  scatter  in  the  derivatives  is  excessive 


Because 

of  the  variance  in  all  derivatives,  especially  the 
angle- of-attack  derivatives,  the  angle-of -attack  derivatives 
were  not  estimated  in  subsequent  runs.  These  derivatives 
are  normally  neglected  when  decoupling  the  modes  of  motion 
into  a longitudinal  set  and  a lateral-directional  set.  The 
side-f orce- surf ace-to-rudder  interconnect  prevented 
estimation  of  the  side-force  surface  and  rudder  control 
derivatives  independently.  A second  test  flight  was 
proposed  using  'no  interconnect  between  the  surfaces.  In  tne 
absence  of  new  data,  the  rudder  control  derivatives  were 
estimated  from  rudder  alone  input  time  histories  and  these 
derivatives  were  then  fixed  at  estimated  values.  The 
side-force  surface  derivatives  were  theD  estimated  from  the 
Flight  1 side-force  surface  time  histories. 

The  two  aileron  input  time  histories  produced 
reasonably  consistent  stability  derivative  estimates.  This 
suggested  that  separating  estimates  by  sets  according  to  the 
type  of  control  would  reduce  the  variance  in  some  estimates. 
The  information  gained  from  each  set  was  incorporated  into 
the  -initial  estimates  for  the  next  set.  Table  5 lists  the 
seguence  of  time  histories  used  and  the  derivatives  which 
were  held  fixed. 

U - 
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TABLL  ,5.  SEQUENCE  OF  TIKE  HISTORIES  USED  IN  ESTIMATING  THE 
LATERAL-DIRECTIONAL  DERIVATIVES. 


CONTROL 

NUMBER  OF  DERIVATIVES 

TIKE  HISTORIES  HELD  FIXED 

AILERON 

7 - 

RUDDER 

“ S <ci  -co 

ea  A6a  6a 

FLT.1  SIDE 
FORCE  SURFACE 

3 CI  'CY  ,Ci  ,C1  ' 

6a  er  J’6a  6r 

C , C 

d n 

6a  er 

FLT.2  SIDE 

6 

FORCE  SURFACE 

The  side-force  coefficients,  C^  and  C were  not 

P Tr 

calculated  because  the  contribution  of  angle-of-attach  to 

the  dimensional  quantities,  Y and  Y , could  not  be 

P r 

accurately  removed. 
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4.2  Estimates  Using  Aileron  Input  Tine  Histories 

Seven  aileron  pulse  tine  histories  with  deflections  of 

7-deg  or  more  were  used  to  estimate  the  VRA  stability  and 

control  derivatives.  Table  4 summarizes  the  average 

estimates  and  associated  standard  deviations  expressed  as 

precentages  of  the  mean  values.  The  variances  of  all 

stability  derivatives  except  C were  reduced.  Only  C.  had 

“r  p 

a standard  deviation  less  than  10  percent  however. 

The  Dutch  roll  mode  was  not  excited  sufficiently  to 
produce  perturbations  in  yaw  rate  and  sideslip  angle  time 
histories  that  had  high  signa 1- to-noi se  ratios.  Maximum 
sideslip  perturbations  were  less  than  3.5  deg  and  maximum 
yaw  rate  perturbations  were  less  than  6.6  deg/sec.  In 
addition,  a dynamic  instrument  error  in  the  sideslip  vane, 
discussed  in  Chapter  3,  makes  side-force  derivative 
estimates  questionable.  .Although  the  estimated  lower  bounds 
of  the  standard  deviations..  (CramSr-Bao  bounds)  of  these 
traditional  Dutch  roll  parameters  are  lower  than  the 
standard  deviations  obtained,  it  must  be  remembered  that  the 
actual  time  histories  and  the  analytical  model  are  subject 
to  measurement  errors,  turbulence  inputs,  and  modeling 
errors.  The  presence  ~of  modeling  errors  is  recognized 

4 - 
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immediately  when  the  estimates  are  grouped  by  direction  of 


aileron  pulse.  The  variances  in  C , C , C and  C 

*r  *0a 

are  reduced  when  estimates  are  separated  into  ' groups  by 
direction  of  aileron  pulse.  These  estimates  are  separated 
by  more  than  two  standard  deviations  (Appendix  A).  A 
5-percent  aileron  scale  factor  correction  that  represents  a 
more  linear  response  of  the  ailerons  did  not  significantly 
reduce  the  separation  of  the  estimates.  tiore  research  is 
needed  to  investgate  this  seemingly  nonlinear  response  to 
aileron  inputs.  A nonlinear  model  would  be  needed  to 
identify  an  asymmetric  aircraft  response. 

Figures  6 and  7 illustrate  some  of  the  difficulties  in 
using  aileron  time  histories.  Time  history  matches  of  roll 
rate  and  bank  angle  generally  are  aood.  Sideslip  angle 
matches  are  poor,  reflecting  the  dynamic  error  in  the 
measurements  and  the  low  siqnal-to-noise  ratio,  and  possibly 
a poor  match  for  the  Dutch  roll  damping  derivative,  Cq 
The  lateral  acceleration  matches  are  fair,  usually  following 
gross  trends  well,  but  the  noise  in  the  control 
measurements,  especially  for  the  side-force  position,  is 
noticeable  in  these  time  histories.  A heading  angle 

measurement  might  be  used  in  place  of  the  sideslip 
measurement . 
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FIGURE  6.  COMPARISON  OF  FLIGHT  1 DATA  ( K ) WITH  TIKE 
HISTORIES  (E)  COMPUTED  USING  THE  PRESENT  STUDY  STABILITY  AND 
CO  NT  SOL  DERIVATIVES  FOR  A POSITIVE  AILFRON  PULSE  INPUT. 
(ABSCISSA  IS  TIME  IN  SEC) 
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FIGURE  7.  COMPARISON  OF  FLIGHT  1 PAT A (M)  KITH  TIME 
HISTORIES (E)  COMPUTED  USING  THE  PRESENT  STUDY  STABILITY  AND 
CONTROL  DERIVATIVES  FOR  A NEGATIVE  AILERON  PULSE  INPUT. 
(ABSCISSA  IS  TIME  IN  SLC) 
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It  was  expected  that  rudder  inputs  would  better  excite 


the  Dutch  roll  mode  and  thus  reduce  the  variance  in  several 
stability  derivatives. 

4. 3 Estimates  Using  Fudder  Input  Time  Histories 

Four  rudder  doublet  (+ 2 . 8-deg  to  ♦ 6 . 3-dey ) input  time 
histories  produced  stability  and  control  derivatives  with 
reduced  variances  in  most  cases.  Tanle  4 summarizes  the 
estimates  using  these  time  histories.  The  major  reductions 
in  variances  were  in  derivatives  whose  associated  state 
perturbations  were  much  larqer  than  those  associated  with 

the  aileron  input  time  histories  (C  , C.  , C ) . 

1 fi  r n r 

In  order  to  determine  the  best  use  of  the  estimates 
obtained  from  the  aileron  input  time  histories,  three  runs 
were  made  evaluatinq  the  benefit  of  weighting  changes  in 
estimates  from  a priori  values.  The  first  run  used  as 
starting  values  the  estimates  obtained  from  the  aileron 
input  time  histories.  Ihe  second  run  also  used  these 
startinq  values  and,  in  addition,  weighted  excursions  from 
these  values  using  the  following  equation  to  calculate  the 
weignting  term  for  the  D2  matrix  (Eg  20): 
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J 


d 

n - 

i JL1L 

i j 

2 

(40) 

D 

a 

i j 

2 

where  c 

is 

the  variance  in 

the 

ith 

row/jth  column 

i j 

der i va  ti ve 

of 

the  A or  B matrix. 

an  d n 

is 

the  number  of 

aileron 

input  time  histories 

used 

i n 

estimating  that 

derivative.  A third  run  used  the  same  starting  values  and  ^ 
priori  weighting,  holding  the  aileron  control  derivatives 
fixed  at  the  previous  estimates  (i.e.,  the  average  or  the 
estimates  for  both  directions).  Since  these  estimates  were 
ontainea  using  aileron  inputs  with  much  superior 
siq na 1- to- noise  ratios  by  comparison  to  the  deflections 
obtained  during  the  rudder  input  time  histories,  these 
estimates  were  the  best  obtainable  set.  Tnis  conclusion  was 
confirmed  by  the  large  departure  of  the  aileron  control 
derivatives  from  their  previously  estimated  values  tor  the 
first  two  runs,  and  also  the  four-  to  five-fold  increase  in 
the  estimated  lower  bounds  of  the  standard  deviations  of 
these  derivatives  (Cramer-Kao  bounds). 

Side-force  derivatives  for  Runs  2 and  3 changed  less 
than  ten  percent  from  estimates  obtained  from  Run  1. 
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I 


Rolling  moment  derivatives  changed  from  30  to  50  percent  for 

Fun  2 and  from  50  to  250  percent  for  Run  3.  All  changes 

were  towards  the  previous  estimates,  and  the  largest  changes 

were  for  L as  expected.  Cram^r-Rao  bounds  increased 
6a 

somewhat  for  these  derivatives,  but  it  must  be  remembered 

that  tnese  bounds  represent  the  lower  bound  of  the  standard 

deviations  obtainable  from  the  time  history  used.  The 

previous  estimates  came  from  time  histories  vitn  lower 

Cram€?r-Rao  bounds  than  those  for  Run  1,  so  these  estimates 

have  higher  confidence  levels. 

The  yawing  moment  derivatives  for  Run  2 were  within  ten 

percent  of  those  for  Run  1.  N For  Rur.  3,  decreased  70 

P 

percent  towards  its  initial  estimate.  N increased  51 

P 

percent  and  K decreased  31  percent.  N increased  only  5 
r er 

percent.  For  all  four  derivatives,  however  the  Cramer-Rao 

bounds  decreased  from  Run  1 to  Run  3 (11  percent  decrease 

for  K ),  indicating  that  the  confidence  in  the  derivatives 
6r 

increased  by  fixing  the  aileron  control  uenvatives. 

Estimates  were  separated  into  groups  having  tne  same 
initial  rudder  deflection.  Reductions  in  variance  and 
estimates  separated  oy  two  standard  deviations  were  obtained 
for  C ' • ci  • cn  aDd  cn  (Appendix  A).  These 

results  again  indicate  that  the  VBA  roll  response  is 
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nonlinear  and  that  its  yawing  response  may  be  as  well.  Bore 
rudder  input  data  would  aid  in  establishing  a yawing 
nonlinearity.  These  results  do  reenforce  the  need  for 
research  on  an  improved  analytical  model. 

There  was  little  or  no  degradation  in  the  time  history 
matches  from  Run  1 to  Run  3.  All  state  time  histories  were 
matched  very  well,  and  the  sideslip  angle  and  lateral 
acceleration  matches  were  far  superior  to  the  lower 
amplitude  aileron  input  time  history  matches.  The  dynamic 
error  in  the  sideslip  angle  measurement  was  not  apparant  for 
angle-of-attack  perturba tionsi of  4 degrees  and  pitch  rate 
perturbations  of  3.5  degrees/second. 

Figure  8 shows  matches  for  a rudder  input  time  history 
using  the  final  set  of  estimates  from  this  study.  Because 
sideslip  angle,  yaw  rate,  and  lateral  acceleration 
perturbations  are  much  larger  than  those  for  aileron  time 
histories  (7-10-fold  increase)  these  matcnes  are  much  better 
than  the  previous  matches.  With  the  Dutch  roll  mode-related 
stability  derivatives  and  rudder  control  derivatives  better 
identified,  the  side-force  surface  input  time  histories  were 
used  for  parameter  estimation. 
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FIG03E  8.  (CONTINUED). 


i*  -20 


I 


6.569  ■» 

1 

1 

rm 

i 

Ml  M 

--■* 

1 

1 

1 

1 

1 

1 

+ 

1 

1 

1 

1 

4- 

1 

1 

1 1 

1 1 

1 

1 

R 

1 

1 

1 

1 

1 

1 

1 

1 

1 1 

1 | 

1 

U 

1 

1 M 

1 

1 

1 

1 

1 

1 | 

D 

1 

1 

1 

1 

1 

1 ‘ 

1 

1 

1 | 

D 

1 

1 

1 

1 

1 

1 

1 

1 

1 | 

E 

1 

1 1 

1 

nrrm 

1 

1 

1 

1 

1 | 

R 

1 

1 

in 

1 

1 

1 

1 

1 1 

| 

1 

1 

1 

i i 

1 

1 

1 

1 1 

ft 

1 

1 

1 

i 

1 

rmrti 

1 

1 

1 1 

1 

N 

0.332  ^ 

--H M — 

--4— 

— 4 r 

-4-M 

4 

4. 

4 

G 

im  n 

1 

1 

i n 

1 M 

1 

1 

irm 

1 1 

m 

L 

l rttn 

1 

1 

i 

1 

1 M 

l rmm  m mi  rm 

uttri  i 

E 

l 

1 

1 M 

i 

1 M 

i rm 

I ri 

i 

rm  rrmtrr  im 

M 1 

l 

1 

1 

1 M 

im 

I 

tmttti 

i 

i l 

i 

D 

l 

1 

1 

1 M 

i 

l 

1 

i 

l l 

i 

E 

1 

1 M 

1 

1 M 

rm 

i 

l 

i 

l l 

i 

G 

l 

1 

1 

i m i 

I 

l 

i 

l I 

i 

1 

M 

1 M 

1 

i 

I 

l 

l 

i l 

i 

l 

1 

mi  m 

1 

l 

1 

l 

i 

l i 

0.0 

2^00 

4.000 

6.000 

0.000 

10.0 

1.398  -l 

+ 

+— 

4. 

— — 4- 

4--M 

— 4— 

M 

M 

ri 

m n 


-0.227 

0 

193.353 


M M 

m rm  n 


•* M-M-  + + M M— r— 


M 

n m 


ri 

m ti 
M M 


m m 
n 

M i 

m n 

M M 


M 

n n 


n 

n m 


ri  m 
im  rrm  n 

185.704  4-M ' 


m n 
n n 


n 

m 

M 

:m 

— +-m-- 
2.000 


n 

4 ri- 

fi 


n — ■*— 


n rtti 
n rm 


FIGURE  H.  (CONCLUDED) 
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4 • 44  Estimates  Using  Side-Force  Surface  Input  Tice  Histories 

Nine  side-force  surface  input  time  histories  were  used 
to  estimate  the  VRA'S  stability  and  control  derivatives. 
Three  were  previously  used  during  the  initial  runs.  These 
were  from  Flight  1,  in  which  a side-force-surf ace-to-r udder 
interconnect  was  used.  Six  were  from  Flignt  2,  for  which 
the  interconnect  was  removed.  No  measurement  of  roll  rate 
was  available  on  Flight  2 due  to  an  instrument  failure,  and 
the  lateral  acceleration  measurement  was  poor. 
Conseguentl  y , the  estimates  were  degraded.  Table  <4 
summarizes  the  results  of  these  estimates. 

Aileron  and  rudder  control  derivatives  were  held  fixed 
for  the  Flight  1 time  histories.  starting  values  for  the 
derivatives  to  be  estimated  were  based  on  estimates  obtained 
from:the  aileron  and  rudder  inputs  time  histories.  i priori 
weighting  was  used  again,  with  the  weighting  matrix  elements 
being  calculated  in  the  same  manner  as  for  estimates  using 
rudder  input  time  histories. 

The  ^ priori  weighting  for  Flight  2 data  was  based  on 
the  Flight  1 side-force  surface  input  results  for  side-force 
derivatives  as  well  as  previous  estimates.  All  control 
derivatives  were  estimated. 
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estimate 


For  these  two  sets  of  time  histories  the 

'p 

is  poor  (Table  4)  due  to  low  sig.na  1-to-rioise  ratios  for 

sideslip  angle  and  lateral  acceleration.  The  largest 

sideslip  angle  perturbation  for  any  of  the  nine  time 

histories  was  7.5-deg.  This  compares  with  the  maximums  of 

from  10  to  15  deg  which  were  obtained  with  rudder  doublet 

inputs.  Greater  care  in  input  design  must  be  taken  to 

accrue  the  maximum  benefit  from  the  side-force  surfaces  in 

identifying  side-force  stability  derivatives. 

The  side-f orce-surface-to-aileron  interconnect  used  on 

Flight  1 effectively  cancelled  side-force  surface  rolling 

moments,  since  Flight  1 time  histories  did  not  excite  the 

roll  mode  to  sufficient  amplitudes.  Cramer-Rao  bounds  for 

rolling  moment  stability  derivatives  usually  were  greater 

than  14  percent  of  the  estimates, with  1 nounds  greater  than 

r 

75  percent.  This  indicates  a low  confidence  level  for  these 

estimates.  No  roll  rate  measurements  were  available  from 

Flight  2 time  histories.  Cramer-fiao  bounds  all  were  greater 

than  10  percent  of  the  estimates  of  1 , L and  L . These 

P P r 

derivatives,  therefore,  were  not  determined  from  tne 
side-rorce  surface  input  time  histories  available. 

Yawing  moment  stability  derivative  estimates  were  not 
as  accurate  as  was  expected.  was  identified  reasonably 
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well,  with  standard  deviations  near  the  Cramer-Rao  lower 

bounds.  The  standard  deviation  for  C from  Fliqht  1 time 

°P 

histories  is  surprisingly  small,  in  fact.it  is  smaller  that 
the  estimated  Cramer-Rao  bound.  The  lack  of  roll  rate 
measurements  and  the  poor  lateral  acceleration  measurements 

for  Fliqht  2 contribute  to  the  variance  in  C and  C 

n n 

P r 

The  side-force  surface  control  derivative  estimates 

from  Flight  2 reflect  the  estimates  from  Flight  1 through 

the  D2  weighting  matrix  (Eg  28).  In  addition,  they  are  not 

perfectly  correlated  with  aileron  and  rudder  derivatives. 

The  C estimate  is  therefore  accepted  over  the  estimate 

1 esf 


from  Fliqht  1 time  histories.  The  very  large  variance  in 


' 6 S f 


the  corresponding  large  variance  predicted 


by 


Cramer-Rao  bounds,  and  the  near  zero  estimate  indicate  that 
this  derivative  was  not  determinable  from  the  data 
available.  Again,  the  side-f orce-surface-to-aileron 
interconnect  used  during  Flight  1 and  the  failure  to  measure 
roll  rate  during  Flight  2 contrioute  to  tne  large  variance. 

Fiaures  9 and  10  are  Flight  2 time  history  matches 
using  present  study  derivatives.  In  Figure  9,  a b-sec 
side-force  surface  pulse  produces  large  state  perturbations 
and  the  matches  are  reasonably  good.  The  sideslip  angle 
measurement  in  Figure  10  suggests  that  the  measurement  error 
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FIGURE  9.  COMPARISON  OF  FLIGHT  2 DATA(K)  WITH  Tift 
HISTORIES  (E)  COMPUTED  USING  THE  PRESENT  STUDY  STABILITY  AN 
CONTROL  DERIVATIVES  FOR  A SIDE  FORCE  SURFACE  PULSE  INPUT, 
(ABSCISSA  IS  TIME  IN  SEC) 
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is  not  corrected  by  the  addition  of  a second  vane.  The 
dynamic  error  is  clearly  visible  and  in  phase  with  the 
side-force  surface  pulses. 

h . 5 Final  Results 

Table  6 compares  the  final  set  of  estimates  with 
unpublished  i analog  matching  results  for  the  VKA  and 
published  estimates  for  a standard  Navion  (without 
side-force  panels).  They  are  presented  to  give  some 
indication  of  the  variance  in  estimates  among  investigators 
and  to  identify  some  effects  of  the  VF  A modifications. 
Reference  ? results  were  estimated  using  a maximum 
likelihood  estimation  approach  similar  to  this  study  on  the 
VHA  prior  tc  the  modification.  The  analog  matching  results 
are  the  most  current  estimates  available  and  offer 
interesting  comparisons. 

Rudder  control  derivatives  for  the  analog  matching 
results  were  not  estimated,  but  they  represent  analytical 
corrections  to  estimates  made  prior  to  the  rudder 

modification.  C.  and  C compare  very  well , however . 

n 

er  “er 

C estimates  from  Ref  3 are  much  higher  than  estimates 
er 

from  other  sources.  A similar  relationship  holds  for  the 
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Analog 
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Ref.  1 

Ref. 18 
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TABLE  6.  (CONTINUED)  . 
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Ref. 3 Ref . 2 

Ref.  1 

Ref . 18 
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-.0865 

-.0485 
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-.06/ 
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-.1556 
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nsa 
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/r  002  /-001 

t005/ 
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cv 

-.0974 

-.  1068 

-.06/  -.08/ 

/-.  08  /-  . 09 

-.08/ 

/-.  08 

-.o  y 

C|*6Sf 

. 0430 

. 0368 

estimate  from  this  study  and  the  analog  matching  result. 

The  large  variance  in  side-force  surface  control  derivatives 

reflects  the  difficulties  encountered  in  this  study,  and 

further  research  is  indicated. 

This  studv  used  the  opposite  sign  convention  from  other 

studies  for  aileron  deflection.  Aileron  control  derivatives 

are  in  poor  agreement  with  analog  matching  results.  C is 

*6  a 

somewhat  larger  in  magnitude  than  the  analog  result  (17 

percent) . Cn  is  5 times  as  large  as  the  analog  result  and 
ea 

is  opposite  in  sign  (due  to  opposite  sign  convention) . 
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Additional  runs  were  made  using  ^ priori  weighting  on  the 
stability  derivatives,  holding  rudder  and  side-force  surface 


der 

i vat 

ives 
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and  allowing  aileron  control 
i do  'S  priori  weighting.  C 

*6a 


6a 


C estimate.  Reference  3 fixed  both  C. 
n 1 

ea  fia 


and  C 


at  wind 


6a 


tunnel  values  in  order  to  obtain  better  estimates  of 

1. 


parameters  correlated  with  these  control  derivatives  ( C,  , 


C ).  Additional  runs  holding  aileron  control  derivatives 

n 

F 

fixed  at  analog  matching  values  resulted  in  poor  sideslip 
angle,  lateral  acceleration,  and  yaw  rate  time  history 


,and  C estimates  doubled  in 
n n 

P r 


matches.  C.  , , • c 

P r 

magnitude  (Appendix  A)  although  C did  decrease. 


Additional  data  using  different  aileron  inputs  are- needed  to 
improve  the  rolling  moment  derivatives. 


4.6  Conclusions 


Initial  estimates  demonstrated  the  dynamic  error  in 
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sideslip  angle  measurements,  and  the  need  for  large  control 


inputs  and  state  responses.  The  scatter  in  estimates  was 
reduced  when  estimates  were  separated  by  type  and  by  initial 
deflection  of  control  surfaces.  Control  interconnects 
prevented  accurate  side-force  surface  derivative  estimates 
for  Flight  1,  and  lack  of  a roll  rate  measurement  and  poor 
lateral  acceleration  measurements  degraded  these  estimates 
on  Flight  2.  An  improved  analytical  model  and  reduced 
measurement  noise  is  needed  to  investigate  the  nonlinear 
roll  and  yaw  responses  documented  in  this  study.  Reference 
3 and  5 suggest  model  improvements,  while  optimal  filtering 
and/or  smoothing  could  reduce  measurement  noise  in  the 
current  data.  Heading  angle  information  might  be  used  in 
place  of  (or  in  addition  to)  sideslip  angle  measurements  to 
prevent  biases  of  the  side-force  derivative  estimates  due  to 
instrument  errors. 

The  poor  agreement  of  aileron  control  derivatives  with 
analog  matching  results  indicates  the  need  for  further 
research  in  this  area.  Attempts  to  improve  the  aileron 
derivative  estimates  in  the  present  research  included  fixing 


lateral 

poorer 


acceleration 
than  those  for 


time 

history 

other 
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matches  are 
deflections . 


generally 
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5.  LONGITUDINAL  DERIVATIVE  DETER MI NATI ON 

The  VRA  longitudinal  derivatives  were  estimated  using 
eight  elevator  doublet  input  time  histories.  The  measured 
guantities  used  were  a ngle-of-attack , velocity,  pitch  angle, 
pitch  rate,  normal  acceleration,  and  elevator  deflection. 
Parameter  estimates  from  these  time  histories  were  averaged, 
and  their  variances  indicate  the  quality  of  the  estimates. 
Attempts  to  reduce  the  variances  included  fixing  single 
derivatives  from  pairs  of  highly  correlated  derivatives  and 
separating  results  by  initial  direction  of  the  elevator 
deflection.  The  effects  of  poor  signal-to-noise  ratios  and 
nonlinear  aircraft  responses  were  documented  in  this  way. 

Initial  runs  clearly  demonstrated  the  dynamic  velocity 
measurement  error  due  to  the  location  of  the  aircraft  static 
ports.  Attempts  to  estimate  longitudinal  force  derivatives 
resulted  in  minimization  algorithm  divergence  due  to 
inadequate  velocity  perturbations  produced  by  aircraft 
response  to  elevator  doublets.  For  these  reasons,  the 
velocity  equation  in  Eq  7 was  neglected,  and  the  associated 
stability  and  control  derivatives  were  not  estimated. 


5 


1 


5.1  Initial  Puns 


Seven  time  histories  (each  of  5 sec  duration)  were  used 
to  estimate  the  stability  and  control  derivatives  associated 
with  the  short  period  mode.  Table  7 summarizes  the  results 
from  these  initial  runs  along  with  standard  deviations  as  a 
percentage  of  the  estimates.  Two  immediate  difficulties  are 
evident  in  these  results.  The  first  is  that  the  sigD  of  C 

It 

is  inconsistent  with  theory  and  its  standard  deviation  i 
large.  The  second  is  that  the  relative  magnitudes  of  C 

7ee 

and  C are  inconsistent  with  theory. 

'"fie 

TABLE  7.  LONGITUDINAL  DERIVATIVE  ESTIMATES  AND  PERCENT 
STANDARD  DEVIATIONS  FROM  ELEVATOR  DOUBLET  INPUTS. 


Initial 

Runs 

Final  Runs 

c 

ZQ 

-3. 6565 

7.5 

-5.6226  12.7 

cz 

q 

18. 9338 

86.5 

Cz6e 

-1 . 1496 

27.5 

C 

a 

a 

-.4068 

24.9 

-.6167  19.8 

C 

"q 

-18.5956 

32.7 

-7.8665  46.9 

c 

•fie 

-1.2776 

25.5 

-.8207  21.0 
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cr  t/> 


Reference  3 suggests  that  C is  an  important  parameter 

q 

for  tne  Navion.  Since  the  addition. of  side-force  panels  and 

the  increase  in  rudder  area  should  have  little  effect  on  the 

longitudinal  stability  derivatives,  C was  included  in  the 

L 

<3 

analytical  model  and  estimated  to  be  positive.  A positive 
C is  not  consistent  with  aerodynamic  theory.  The 

zq 

incremental  normal  force  produced  by  positive  pitch  rate 
perturbations  is  primarily  due  to  an  increase  in  the  angle 
of  attack  at  the  horizontal  tail.  An  increase  in  angle  of 
attack  results  in  a negative  normal  force.  The  positive  C 

L* 

q 

estimate  is  also  inconsistant  with  the  negative  C estimate 

’ m 

q 

(Ref  5)  . The  large  variance  in  the  estimates  of  this 
derivative  indicates  that  it  is  not  readily  determinable 
from  the  available  data.  Different  control  inputs,  e.  g., 
direct  lift  pulses  from  wing  flaps,  might  help  isolate  this 
derivative. 

The  pitching  moment  due  to  elevator  deflection  is 
primarily  due  tc  the  change  in  normal  force  of  the 
horizontal  tail  due  to  elevator  deflection: 


5 


3 


the  V?A;  however,  the  present  parameter  estimates  do  not 

conform  to  this  theoretical  relationship. 

The  normal  acceleration  time  histories  from  these  runs 

initially  show  the  effect  of  the  erroneous  C and  C 

Zg  Zfie 

The  estimated  time  history  initial  response  to  the  elevator 
deflection  is  a visible  acceleration  opposite  to  the 
intended  acceleration  as  seen  in  Figure  11.  A pitch-up 
elevator  deflection  initially  prgduced  a large  upward 
acceleration  of  the  c.g.  in  the  estimated  time  history  (I). 
This  response  was  clearly  not  evident  in  the  measured  normal 
acceleration  (*) . 

Reference  3 jand  10  suggest  that  the  inverse  of  the 
second  gradient  of  th£:cost  function  with  respect  to  the 
unknown  parameters,  as  calculated  by  the  minimization 
algorithm,  is  a useful  approximation  to  the  parameter 
covariance  matrix.  Table  8 shows  this  matrix,  normalized  by 
the  diagonal  elements.  The  off— diagonal  terms  of  this 
matrix  are  approximations  to  the  correlation  coefficients  of 
corresponding  parameters.  The  high  correlation  between  C 

Z6e 

and  is  expected  from  their  geometric  relationship.  In 

69 

addition,  the  aerodynamic  bias  estimates  are  highly 
correlated  with  the  control  derivatives  possibly  because  of 
a bias  in  the  trim  elevator  deflection  measurement. 
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(E)  CONFUTED  USING  INITIAL  ESTIMATES  OF  THE 

AND  CONTROL"  DERIVATIVES  FOR  AN  ELEVATOR  INPUT. 
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FIGU-r  11.  (CONTINUED). 
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TABLE  8.  ESTIMATED  LONGITUDINAL  PARAMETER  COVARIANCE  MATRIX 
NORMALIZED  BY  THE  DIAGONAL  ELEMENTS. 


■ 


not  determinable  from  the  data.  In  addition,  the 

aerodynamic  biases  were  calculated  to  balance  the  equations 

of  motion  initially  and  were  held  fixed  during  the 

estimation.  The  results  of  this  second  set  of  runs  are 

summarized  in  Table  7.  Note  that  the  C estimate  increased 

Z 

a 

in  magnitude.  The  short- period  mode  damping  derivative,  C 

B 

g 

decreased  in  magnitude.  c is  now  approximately  three 

*6e 

times  the  magnitude  of  C 

ee 

In  attempting  to  reduce  their  standard  deviations,  the 
estimates  were  separated  into  sets  by  direction  of  initial 


elevator  deflection.  The  variance  in  C and  C decreased 

m m 

« g 

for  the  estimates  separated  in  this  way  and  all  estimated 
derivatives  had  significantly  different  mean  values  based  on 
this  separation  (see  Appendix  B) . This  nonlinear  response 


to  elevator  inputs  again  suggests  that  the  analytical  model 
should  be  reexamined  and  revised.  It  also  suggests  that 
smaller  amplitude  inputs  must  be  used  to  remain  witnin  the 
linear  range  of  aircraft  response.  The  penalty  of  poorer 
signal-to-noise  ratios  on  the  guality  of  estimates  dictates 
that  the  measurement  noise  must  be  reduced  to  maintain  high 
signal-to-noise  ratios. 

Dse  of  a Kalman  filter/smoother  during  data  processing 


is  one  approach  to  reducing  noise  in  the  data  provided  by 


5 


8 


the  VRA's  existing  sensors.  Care  Bust  be  exercised  in  this 
procedure  to  prevent  parameter  estimates  from  being  biased 
by  the  dynamic  model  assumed  for  the  Kalman  filter/smoother. 
A long-term  approach  to  noise  reduction  is  the  replacement 
of  the  hardware  used  for  data  processing  in  this  study  with 
equipment  better  suited  for  digital  data  analysis. 

Table  9 compares  the  estimates  derived  in  this  study 
with  analog  matching  results  and  published  estimates  for  the 


Navion  in  different  configurations. 


is  higher  than 


expected  from  earlier  results.  Reference  3 suggests  that 
failure  to  include  C in  the  analytical  model  results  in  a 

Zq 

larger  amplitude  C estimate.  The  results  from  the  initial 

Za 

runs  suggest  that  a positive  C is  required  to  reduce  C 

Lj  _ / 

g a 

estimates  from  the  given  data.  Additional  runs  including 

C ^ in  the  model  resulted  in  increased  magnitude  C z 

q a 

estimates  for  negative  C estimates.  The  variance  in  the 

zq 

C estimates  still  was  very  high.  Additional  data  would 

q 

better  determine  the  relationship  of  C and  C . 

z z 

q « 

The  three  estimated  aoaent  derivatives  have  aagnitudes 
smaller  than  the  analog  matching  results,  however,  they  are 
within  one  standard  deviation  of  those  results,  with  the 
exception  of  C , which  is  half  the  analog  value.  The  C 


m 


6e 


estimate  is  approximately  three  times  the  C estimate. 

2ee 
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TABLE  9.  COMPARISON  OF  LONGITUDINAL  DERIVATIVE  ESTIMATES 


WITH 

ANALOG 

HATCHING  AND 

PUBLISHED 

RESULTS. 

Ana  log 
Match 
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Ref  .2 

Ref . 1 Ref. 18 
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/-I.  5 
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-1.5 

Additional  data  using  improved  elevator  input  shapes 
and  flap  deflections  could  reduce  the  correlation  of  several 
of  the  derivatives.  Ar  improved  estimation  method  also  is 
needed  to  reduce  the  biases  in  estimates  resulting  from 
holding  elements  of  correlated  sets  of  derivatives  fixed. 

Time  history  matches  using  the  stability  and  control 
derivatives  estimated  in  this  study  show  good  agreement.  In 
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Figure  12,  the  angle-of-attack  matching  and  the  normal 
acceleration  matching  do  suggest  some  inconsistency  in  the 
measurements.  while  the  estimated  time  history  tends  to 
undershoot  extrema  in  the  angle-of-attack  measurements,  it 
tends  to  overshoot  extrema  in  the  normal  acceleration 
measurements.  This  demonstrates  the  effectiveness  of  the 
minimization  algorithm,  as  well  as  the  adverse  effects  of 
scale  factor  and  bias  errors  in  the  measurements.  An  error 
in  the  upwash  correction  factor  of  the  angle-of-attack 
measurements  might  result  in  the  inconsistency  between 
measurements  sources. 

It  was  suggested  that  an  error  in  the  upwash  correction 
for  anqle-of-attack  measurements  also  might  result  in  a C 

Zq 

estimate  with  large  variance.  Six  additional  runs  were  made 
using  two  elevator  doublet  input  time  histories  (Appendix 
B) . These  runs  included  in  the  model  and  varied  the 

q 

upwash  correction  factor  from  1 (no  correction)  to  1.9.  The 
results  suggest  that  there  is  a strong  relationship  between 
the  upwash  correction  and  C but  that  the  variance  in  the 

zq 

C estimate  is  not  reduced,  as  seen  in  Figure  13.  While 

zq 

some  matches  improved,  the  undershoot-overshoot  relationship 
was  never  reversed  ( a measurements  were  consistently 
undershot  by  a estimates)  . 


f 
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FIGURE  13.  THE  RELATIONSHIP  BETHEEN  ANGLE-OF-ATTACK  OPbASH 
CORRECTION  FACTOE  AND  THE  COEFFICIENT  FOR  THE  RATE  OF  CHANGE 
OF  NORHAL  FORCE  WITH  RESPECT  TO  PITCH  RATE  FOR  TWO  FLIGHT  1 
ELEVATOR  INPUT  TIHE  HISTORIES. 


Aerodynamic  biases  were  held  fixed  during  these  runs. 

Tne  four  data  points  at  a correction  factor  of  1.6  show  that 

failure  to  use  the  correct  analytical  model  results  in 

biased  estimates.  Two  of  the  data  points  (triangles)  are 

the  estimates  from  initial  runs  for  these  time  histories. 

Here  Z , Z , and  the  aerodynamic  biases  were  estimated  as 
q ee 

well  as  the  other  short  period  mode  derivatives.  The  other 
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two  points  (squares)  were  estimated  with  the  aerodynamic 

biases  and  Z being  fixed.  This  result  clearly  indicates 
ee 

the  need  for  further  research  on  the  analytical  model. 

Additional  runs  were  made  to  determine  the  effects  of  a 
non-zero  pitching  moment  with  respect  to  rate  of  change  of 
angle  of  attack  (C  ) . This  term  was  introduced  using  the 

I • 

a 

acceleration  transformation  matrix,  fi: 


Rx  = A x + B u (42) 


Inverting  B and  premultiplying  the  A and  B matrices 

correctly  accounts  for  the  unsteady  aerodynamic  effect  on 

the  stability  and  control  derivatives  of  this  term.  CB  was 

a 

set  at  -4.3e5,  the  value  from  analog  matching  results.  This 

C increased  C +C*  by  25  percent  and  increased  C by 

a q : a 6e 

ID  percent  towards  the  analog  matching  results.  A similar 

result  is  reported  "in  Ref  3.  changed  by  less  than  2 

a 

percent  froir  its  previously  estimated  value. 


5 . 3 Conclusions 


Initial  estimates  demonstrated  the  dynamic  velocity  error 
and  the  small  velocity  perturbation  response  to  elevator 
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doublet  input  time  histories  of  5-sec.  In  addition, 
physically  nonrealizable  estimates  of  certain  derivatives 
and  prediction  of  high  correlations  in  pairs  of"  estimates 
suggested  that  the  associated  derivatives  should  be  held 
fixed  during  estimation.  These  runs  produced  good  time 
history  matches  and  fair  comparisons  with  previous  work. 
Separation  of  results  by  initial  deflection  reduced  the 
variance  in  some  estimates  and  clearly  indicated  a nonlinear 
aircraft  response.  This  recommends  further  refinement  of 
the  analytical  model  to  include  nonlinear  terms  and  use  of 
lower  amplitude  elevator  inputs  to  stay  within  the  linear 
range  of  the  lift  curve  slope.:.  This  in  turn  requires 
reducing  measurement  noise  eitner  by  smoothing  existing  data 
or  by  employing  improved  sensors.  The  effects  of  scale 
factor  errors  in  angle-of-attack  measurements  were 
investigated,  and  the  results  indicate  that  further 
verification  of  the  upwash  correction  should  be  made.  In 
addition,  an  improved  method  for  estimating  highly 
correlated  derivatives  must  be  devised.  Inclusion  of  a 

fixed  value  of  CM  seemed  to  improve  moment  derivative 

a 

estimates  and  should  be  estimated.  This  again  would  require 
modification  of  the  analytical  model. 
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6.  CONCLUSIONS  AND  RECOBMENDATIONS 

The  research  for  this  thesis  resulted  in  several 
conclusions  and  recommendations  about  analytical  models, 
data  collection  and  derivative  estimation: 

. The  absense  of  measurement  and  process  noise  in  the 
flight  data  is  not  always  a reasonable  assumption,  as 
the  quality  of  estimates  is  directly  related  to  the 
validity  of  this  assumption. 

. Instrument  errors  result  in  the  degradation  of  the 
quality  of  derivative  estimates.  The  aircraft  static 
ports,  angle-of-attack  and  sideslip  vanes,  and  lateral 
accelerometer  should  be  studied  to  determine  the 
effects  of  high  frequency  inputs,  rapid  changes  in  air 
flow,  and  upwasb  effects. 

. Calibration  procedures  should  be  reexamined  with 
the  aim  of  reducing  human  error.  Simple  calibration 
instruments  can  be  constructed  that  result  in 
consistent  results. 

. Improved  control  input  designs  can  now  be  verified 
using  the  VRA  digital  flight  control  system,  which  is 
capable  of  generating  repeatable  inputs  of  arbitrary 
pulse  (or  wave)  shape. 
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. The  replacement  of  analog-to-digita 1 conversion 

i 

hardware  with  eguipment  better  suited  to  digital  data 
analysis  will  reduce  measurement  noise  and  biases  added 
by  the  multitude  of  components  now  used. 

. Mass  and  inertia  information  must  be  determined  for 
the  VRA  in  its  current  configuration.  The  effects  of 
inertia  coupling  should  be  documented  prior  to 
identifying  aerodynamic  coupling  parameters. 

. Poor  signal-to-noise  ratios  aggravated  by  failure 
to  sufficiently  excite  all  modes  of  motion  result  in 
poor  quality  derivative  estimates.  These  effects  are 
visible  in  the  time  history  matches  produced. 

. The  nonlinear  response  in  roll,  pitch, and  yaw  to 
control  estimates  is  documented  by  the  reduction  in 
variance  of  estimates  separated  by  initial  deflection 
of  the  control  surface.  This  suggests  that  the 
assumptions  used  in  deriving  the  analytical  model  must 
be  reexamined  and  additional  effects  must  be  modelled 
in  the  equations. 

. Cramer-Fao  lower  bounds  to  the  standard  deviations 
in  estimates  are  reasonable  confidence  level 
indica tors, a lthough  they  rarely  reflect  true  standard 
deviations. 
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. An  improved  method  for  estimating  highly  correlated 
derivatives  must  be  devised.  Improved  control  input 
desiqn  is  a fruitful  direction  for  research. 

. Nonphysically  realizable  estimates  of  derivatives 
can  cause  obvious  errors  in  time  history  matches  that 
are  not  corrected  by  the  minimization  algorithm. 
Measurement  and  process  noise  can  account  for  these 
estimates  as  well  as  poor  excitation  of  modes  of 
motion . 

. Good  time  history  matches  do  not  guarantee  good 
parameter  estimates  in  the  presence  of  measurement  and 
process  noise.  Hodeling  errors  and  noise  and  biases  in 
the  instruments  degrade  the  quality  of  estimates. 

Further  research  is  needed  on  the  aircraft  analytical 
model  and  minimization  algorithm.  Instrument  studies  should 
precede  further  -estimation  research.  Input  design  studies 
orrer  improved  quality  of  estimates  even  with  current 
measurement  and  process  noise  problems. 
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APPENDIX  A 


i 


TABULATED 


L ATE  PAL -DI RECTI  ON  A L 
ESTIMATES 


DERIVATIVE 


The  Modified  Maximum  Likelihood  Estimation  program  used 
in  this  study  produces  estimates  of  stability  and  control 
derivatives  and  lover  bounds  to  the  standard  deviations  of 
the  estimates  (Cram^r-Rao  bounds).  Table  10  lists  the 
lateral-directional  derivative  estimates-  and  Cram^r-Rao 
bounds  for  the  twenty  time  histories  discussed  in  Chapter  4. 
Zero  values  for  the  Cram&r-Rao  bounds  indicate  that  these 
derivatives  were  held  fixed  at  an  assigned  value.  Mean 
values  and  standard  deviations  are  also  listed. 

Runs  1 through  7 in  Table  10,  are  estimates  from 
aileron  pulse  input  time  histories  for  Flight  1.  Taken 
together,  the  standard  deviations  of  all  stability 
derivatives  and  aileron  control  derivatives  are  reduced  from 
values  for  all  estimates.  Standard  deviations  are  reduced 
further,  for  several  derivatives  by  grouping  estimates  into 
sets  according  to  direction  of  aileron  deflection  as  in 
Table  11  (3  positive  and  4 negative  inputs).  This  reduction 
in  variance  indicates  a nonlinear  response  of  the  VRA  to 
aileron  inputs. 


A 


1 


^j  — — t fsj  in  in  >d  in  o in 1 - c j cti  cti  in  m t 

■3-M3I.  in  in  r-j  — CTi<d  rsi  — ru  — o o o 

in  C'  in  ro  co  — a i "7  in  co  in  in  id  c>  o id  in  r j m 

® in  ci  i . i'_  a1  in  o o ct»  id  cm  — — cm  ct>  ij  co  in 

ruo--noQ-"D“00ooot;'OCi- 

K'OQOOE'OOOOOOOOCiOOOOQ 
+ + + 4+  + + + + + + + T + + + + + + + 
iDioixiD',JorxiDoinin  — rMin<xi^irxixiT  — 

cn  rM  co  co  m co  cn  in  ix  tm  co  co  in  cd  c?  co  lx  ix  — in 

® o — cti  cr>  ix  in  ai  co  ■c  cti  — 1 cm  in  ix  — 1 csj  — ro  co 

ix  ix  in  u ix  ix  c.  in  id  rn  id  in  co  cj  Ci  (\i  an  in  ® 'J 

co  cd  ij  in  cm  — ■7  in  o rx  in  in  in  cj  id  cj  in  co  ix  in 



f . i r~ \ i rr i Pwi  r\ i i*r*  rri  rri  ^i  rr i q i ^ i rr * o Cd  CD  ^ ■ * o O 

l l I I l I l l I 

ro  g . cm  in  — co  in  id  ix  in  cd  cd  o U cm  oc  "7  in  in 
in  ■cr  cd  in  uo  *7  v id  ^7  cd  cm  Ci  o cd  ix  in  co  cn  ^7  ^7 

m co  cm  — c>  — in  in  rn  ix  cd  o C'  o — in  ix  co  cti  cd 

ix  a.  iJ  "7  in  co  cn  — r\j  — i o c d o co  ^7  in  — ro  in 

CM  G'  — ^7  — G'  Ci  Gi  Gi  Gi  Ci  G'  Ci  CD  tx  Ci  ID  Ix  CTI  fx 

Cl  Cl  O CD  C'  C*  C'  C'  CD  C»  Ci  C O C Ci  — Gi  C C rM 

+ + + + + + + + X + + + + + + + + + + + 

— ID  lx  — |x  Ix  CD  '?  — CD  CM  in  CD  — in  — CT  'JCOIx 
ci  in  — — U — in  inirjucocoiDOix«  ix  c ix 
cn  co  in  cd  "7  cn  in  co  id  "3  7 co  — o in  co  — do  do  cm  co 

a id  id  id  7 or  in  ng  ix  in  o — co  cti  cti  cn  co  ix  — cn  — 

0 cd  — v cm  — in  in  in  'J  in  in  in  in  n cn  cn  n in  n 

f\ i iV i Cxi  CD  i-i  fxi  Ci  Ci  G)  i x i rx i r»i  t xi  fTi  Q Cj  O CD  O CD 

1 I I I I 

fM  N-  G>  G ro  <X>  G G G>  G«  G Gl  O O — CO  CD  CD  G 

«—  ro  ro  n-  in  ro  rv  c q d k q c c co  ^co^k  ro 

G CO  CTl  CD  K-  N-  Cl  CD  C CO  O Cl  C CO  C-  fO  — CD  fO 

JO  CD  — • <M  Gi  CD  — • C'  Ci  Ci  Ci  G*  CD  — • f'w  If)  f'-  Ci  CV 

G O G G Q G*  d C.  Q O - O d G ^ I'-  <M  <M  M C 

C*  G*  CD  Ci  Ci  G1  CD  d CD  G CD  CD  d CD  CD  G G'  G*  G* 

+ + + + + + + + + + + + + + + + + + + + 

o ^ M co  ci  — • ro  in  — in  co  c in  o — ro  ro  o — ro 

in  id  — • fv-  cn  co  co  in  m g»  — n_  c-  <m  g cn  cn  co  ro  cn 

gc  co  — • c c»  in  co  ^ 'si  ^ ^ in  ^ in  in  co  in  co  *n  ro 

d»  ro  (MGrOGiGMM-*  — --G  — — — 'CTW  M-  CM  C\J  VO 
Cl  Cl  C Ci  d d d CD  Q d G G O C Q0G^  G»  CM 

CD  G'  G*  CD  CD  G'  G'  CD  CD  CD  CD  C»  CD  G'  CD  CD  ^ CD  CD  CD 

I I . __  I I I 

Cl  CD  CD  Cl  Gi  C'  Ci  G'  CD  G'  CD  G d Q d d Q d d CD 

d G'  CD  G*  Gi  CD  Gi  d CD  i~i  i Gi  CT*  CD  d d CD  "•  i i‘  1 1~ . i i\  i 

G'  G'  G*  G'  G'  G»  C«  G' d d Q Q D D G' d C CD  Q C 

G'  CD  Gi  G'  G'  CD  CD  Q G'  G'  G' d G'  G1  O Q G Q Q Q 

QGOOQOddQGidQQOQGGiCiOCD 

000OG'QQGQ0QQGiQQOGiG)QO 
+ + + + + + + + + + + + + + + + + + + + 
0G0GG00G0G0Q000G0Q0G 
, G CD  G*  CD  CD  CD  G*  CD  Gi  CD  Q CD  G1  G1  G*  CD  C*  Gi  *2?  CD 

G G'  G'  C'  G G'  G G G G G G Gi  Gi  G G Gi  G'  G G 

• Gi  Gi  CD  CD  ^D  CD  CD  d CD  G1  CD  d CD  d CD  d CD  C*  CD  CD 

CHEHSOOQODEIQOQQSiOOQOClOO 

G G G Q.Q  GGGGGGGGG'GGGGGG 

(MMXi(TUnimii0000000  000Q0O 
■ ' Gi  Gi  G'  Gi  G'  Gi  CD  Gi  G*  G*  CD  d C*  C*  CD  CD  CD  CD  CD 
QGGG'GGGGGGGGGGGGGGGG 

Gi  G*  Cd  CD  d C- 1 i\  i r . i i\ i i\  1 C . i t‘ > 1 C-.i  CD  C v i f 1 1 r _ 1 CD  d G) 


+ -4-  + + + + + + + + + + + + + + + + + + 
G'GG'GGGGGGGGGGGGGGGGG 
0 0 0 0 0 0 0 0 0 0 0 0 0 0 00  0 0 0 0 

G'  Gl  Gi  d CD  CD  d d d CD  Gi  i~P  Cxi  i" > i d Cxi  i' . i Cxi  ^D  i'i i 

GC'GC'GGGGGGGGGGGGGGGGG 

GGGGGGGGGGGGGGGGGGGG 

viMioixcoiM  — co  ix  ■'T  in  g>  — id  — miniMcnm 
rniD<M(MiDi0GirDnjcnmcicninEMX>OQcnin 

(m  — 'j-CT-CTco^J(MCT'innjinco  — gj  — — omco 
crnnixcoin^rrx  — — . — pgvmruai  — cooixid 
) — ogiCigpo  ©gigjgigioj  — — om  — — — tm 

i 

00  0 0000000100)0001000  0 0 
+ + + -r  + + + + + + + + + + + + + + + + 

^ CM  CT)  in  — • CT>  CM  G G CM  <M  r%-  CM  ^ in  -SJ  CT)  — • G 

«ncnCTiUDGGro<nroirjroG>—*G--C'joroN  — 
G (M  I'-  ro  G — — (M  ^ ro  — ^ CM  — CM  Gi  G G G ro  G 
hino^ininDrxicim'croMnaiOOin'jMj? 
ui  mu?  in  cc'd  cn  d 0 co  k o ro  UT  d N ID  If)  ^ d CTi 
i m 

* dCiCiOdd’HOQd’-'’-,d'H-M^C'  — *-*01 
I I I I I I I I I I I I I I I I I I I I 
| z — * <m  ro  ^ in  in  r^  co  cn  cd  — cm  ro  •sj  in  g r>_  co  cti  g 


G0\‘ 
H CD  Gl 
<L  G • 
G • G 
G 
CD 

<M  . 

+ 

CD 
CD 
z : cd 
O CD 

a o 

LU 

CO 

<L  — X. 

CD  CD  G 
0?  • 
> G«  CM 
LU  ro  ro 
q • 

CD 
Q + 
h-  ^ 
if)  G 

^r 

Q rT 

z cn 
a • 

G' 
z i 

<L  CD 
G (M 
H 


CL  DEP IVhT I 
P'JM  E'ETh 


•£■  a.  © IJ-,  ©.  N.  fO  CO  G 7 r.OtfJM  - 7 G oj  U"> 
fO  till  LO  I'-  C 'X'  'Jj  '£■  — I - Oj  f - 'jj  fO  IO  .0  7 — ~ Cu 

CO  CTi  Cl  — *7  'O  I _ G CO  — 7 IO  oj  OJ  •£•  rg  r.j  m O *7 

— © — in  © O O G G C Z'  • ■ 

O O G O G O G © © O O O G Z»  O G o O C'  Z' 


G O O G C«  C O C«  © o o c*  o o c o o o o o 

+ + + + + + -f-  + + + -f  + + + + +-f  + -t  + 

r-  cr.  © to  in  r-  ro  r-w  i — • c — ‘£»  rg  rg  co  a.  r - ro  co 

CD  IO  ^N.r.  *£•  — u-'  T-  O rg  Ch  G 7 -f  - IO  © rg  rg 

rv  Uj  7 oj  co  un  co  '£•  rg  •£•  m co  o rg  cr»  rg  u~>  m ro  ro 

u_  in  ro  — oj  co  oj  rg  7 cu  co  '£>  — o o — — o C»  10  — 

OlCiD’-QOQO-OO-QOOOQO  G G CD 

Q 

Ci  Ci  CD  C'  CD  C G C>  CD  O C*  CD  CD  O CD  CD  CD  CD  C'  O 
I III 

in 7 fO  Cl  r - CD  CO  — CD  ‘£>  — O Z*  CD  rg  © *7  C'  CO  UO 

oj  in  cr.  aj  oj  co  — — c«  uj  7 o g ci  cc>  7 g — un  © 

707  - d ro  ro  oj  rg  ro  7 © c o ro  r - o o un  7 

OJ  OJ  — 0*1  OJ  — — CD  CD  CD  C*  CD  C«  CD  Uj  OJ  JO  rg  h-  Cl 

Cl  C'  C'  CD  Cl  CD  Cl  'Z'  C'  C'  © G'  'Z*  CD  — • O C'  — ' © ~ 


Ci  O CD  CD  C'  Ci  C«  O Z'  G O G C«  'Z'  G »Z'  C'  C'  C'  C' 

+ + + + + + + + 7 + + + + + + + + + -+  + 

un  oj  'x'  © co  r-  r_  in  o ro  o co  © ro  in  — r-  oj  co  g 

— rg  ■sT  — oj  7 in  cr.  © un  ro  co  — un  cr-  cd  cj  co  — in 

q:  *7  cd  m un  i n o in  r o ©77©  © r o 7 7 d 7 un 

cu  u-*  *7  oj  rg  co  — rg  rg  o j rg  rg  r-J  oj  rg  o j oj  oj  rg  rg  rg 

G G G G t-J  G G G«  G»  G C«  Z‘  G O C»  G G CD  G G 


G C'  G G G G G G O G G G G G O G G G G G 

I I I I I I 1 

in  I'-  1£|  Cu  ‘jj  CO  Cm  G © © 7 C G G a.  — ro  OJ  CTi  — 

© — © in  — IT,  *7  Cl  O Cl  *7  G G G in  CD  Cg  — ~ 7 

7 U~|  *7  r-'i  ro  OJ  ro  G G C 7 Q G G OJ  " IP  7 (Ji  UIi 

© C<  'Z'  G G G G G G G fO  G G G .O  O J fO  '£'  LTj  G 

G G G G G G G G O G G G G G uGGC'G^ 


G O G G G G G G G G G G G G C G G G G G 

+ + -f  + + + + + 4-  + + + + + + + + + + + 

in  di  to  ‘jj  — ct«  i'-  oj  o ro  co  rg  od  co  ro  co  «£«  r j ro  — 

ro  co  cti  ro  ro  cn  — © r^  *o  in  — oj  CTi  i n cm  co  ud  ro  iri 

<X  CTi  © Oj  • CD  Oj  'jj  Oj  CO  *— • C»  ■©  G — 1 1 '£•  »©  — • rg 

CU  © OJ  G <£'  Lj  «Jj  •£>  Oj  CO  CO  ‘£'  CTi  CTi  G f*''  CTi  *7  L"'  — OJ 


G G G G G G G G G G G G G G G G G G G G 

I I I I I I I l I I I I I I I I I I l l 

© uj  •£•  ro  — ro  Cti  d-  © ro  in  oj  rg  in  © ro  rg  in  co  o 

rg  co  in  *7  © rg  *£*  r _ o »z  un  ©•  rg  — -7  a . *7  co  "7  ld 

— — • o r-'j  u'i  'J  — o r - — c>  o tri  t y-.  txi  r-j  co  — 

K)  r j ro  ro  — ■ ro  c-j  — o — c i co  T co  to  in  tv  *t  o r 

G G 'Zi  G G G G G G G C G C G G G G G C ^ 


G G G G © ©i  G G G G C G G G G C >_J  G G C 
+ + + + + 4-T--+  + + + + + -4_4-  + + + ^‘f 
CTi  un  • £'  — un  un  on  u~  — 01  eg  eg  on  co  co  <2\  '£'■  ”7  u-' 

CO  7 1 0*1  rg  «£i  [ - 7 7 co  — © rg  © Cti  CO  CO  CTi  G — •£' 

to  rg  *7  ( - *7  — r - © *7  T '£  '£'  in  fO  U'i  G '£•  .Ji  T U"' 

7 Cl  r^  cr.  '£•  7 in  CO  OJ  CTi  O J © fOJ'j  I - *7  Cti  o j ro  f-j 

01  OJ  G — G G G G © — G G OJ  fO  —*  G G i— 


©.  ©1  ©1  © ©1  ©1  © © G G G G G C G G G G G O 
I I 

nr  r - un  Z"  rg  eg  — «£'  n co  co  l"i  ro  it*  oj  — ro  oj  eg  n 

- — • - — • ' — ■ - — • • — • — - — • u-'  rg  un  © — — — ©•  in  rg  r - g © 

©1  ©»  ©1  © ©1  ©1  © © ©1  © rg  ©1  — • — rg  ©*  — • 10  [-'.  u_i 


©•  © ©•  © O Ci  Ci  G Ci  G G i_J  G C'  • G G G O 


+ 4 + + + -.-+-*-  + + + + + + + + + + + + 

G r - f w c «£'  ro  »£'  — r-o  — Cti  rg  in  r - •£•  *7  u'i  © lh  l"* 

'£•  CT»  CT«  CO  n G ( - IP  7 a 1 — CO  7 7 — '£»  l-j  »Z»  G 7 

7 U”|  IT,  7 rg  un  IT,  rg  rg  — co  r - Cti  co  a,  rg  |£|  IT,  0*1  0 , 

Q_  U'I  '£|  U'I  U'I  U'  U'I  7 •£'  7 Oj  '£i  7 1-1-7  OJ  7 '£•  f - f - 


iTT  i T T iT  i T T i i " T T T i i T 

— * rg  i£t  ro  r - I - '£•  r - CO  CO  '£'  rg  '£(  U'i  — • to  7 7 U"'  r 1 
ro  f . © — r J CT,  rg  '£1  rg  rg  i£«  f . un  ro  7 O '£•  © Oj  oj 

CTi  un  r _ ' n 7 7 CO  © LD  © 7 CT,  rr,  © 7 r g f*  . rg  — 

© © © c»  c © G G G G ro  r ^ >y  in  ~ *0  tm  7 — 

LD  “ "*•  “ • "*•  “ “ "*•  “ • • * • "■  • * • 


•f  + + + + + + T + + + + + + + + + -+  7 + 

rr.  r.j  f _ «— « u"1  ro  — '£•  7 '.0  7 rg  f - 7 to  un  rg  ro  m O 

'£•  ' D I - CO  U'i  rg  i£i  --  © Oj  Oj  CTi  •'  J — ' 'Ji  7 © 7 G 

U'i  —•Of-  u'»  {O  © — < i£i  U'i  «jj  { _ Oj  7 »D  i »D  I © »j  i »'• , <£i 

•7.  lj,  r . N C'  ( - C - 1 - fO  — ro  © ro  — rg  n U-.  rg  rg 


^ T i j i i iT'iTTTT  i 7 i i i iT 

— r j r-j  7 u'  •£•  r •-  co  a.  c«  — rg  ro  7 un  •£•  r - co  •©  c 


f w r - m co  1 co  i . r . ~7  m a.  — r - *7  co  "*  •£•  m»  in  eg 

i w — to  a.  o i . © r j i '.»  a,  |.*,  co  •/.  o eg  in  a,  u*>  ud  © 

in  T l>-  '£'  I - in  '7  *?  '£•  ud  ( j rv  — — r j r j Ci  • - *7  co 

— C«  C>  o C»  O'  C»  O O G»  b"»  © c»  c*  C'  O'  C*  O'  C*  © 

ci'  C1  C*  O C'  c 'll*  C'  c C«  O'  C'  o C C'  C'  O'  o c*  ©• 

C*  O © G*  O © o z*  o o o o o c»  o o c*  o o c» 

+ + -»-  + + + + + + + + + + + + + + + + + 

w co  in  co  eg  ’7  © co  in  r - a i in  iv  »£»  — c :•  »£■  c*  in  © 

— • r - rg  rs-  k ud  co  rv.  eg  a,  co  »£•  •£•  o in  eg  in  in  r . iv 
V£»  f'j  CT « C* J — • <T ( >£i  rr,  | . © — CO  '£'  *7  — CD  CO  'Jj  *7  rij 

ll  — — o — eg  in  in  © — o eg  *7  in  7 ud  eg  in  v *7  *7 

woo  ©•  o c«  o o o o o c»  o o»  c»  O'  © O O Cl  CD 



C O C'  O © C CD  O'  C O C*  CD  C«  C'  O C'  c<  c»  c c» 
II  II 

^ V 7-KC'i.  rv  iv.  •£'  tj  o o o i •-  in  in  cr.  eg  eg 

CM  '-0  — O in  co  © O — '£  '7  O'  O C'  O CO  eg  UD  '£i  '£j 

to  r , r _ i*_  — * u~*  in  — ~ — in  o o o ud  *7  — c*  ud  cd 

*—  ©>©'  — — G*  O C'  C O CD  © O G U"'  eg  — — eg  '£. 

CZi  CD  CD  CD  CD  C<  CD  C*  CD  C»  O CD  C»  C*  O'  C»  0 — 0 O 

C*  CD  C'  O'  o C'  O'  CD  O'  O O'  O'  O*  O'  O O'  O'  O'  O'  O' 

+ + + + + + + + + 4 + 4 + + + + + + + + 

cr  i in  i - co  'a*  «£•  i - e i in  in  in  *7  — © •£■  u.  lt*  cc»  — r - 

co  e*  j in  co  © l/'»  •£■  iji<  cr.  o eg  ud  © co  ct,  ud  o<  c«  •£•  in 

a:  — 'jd  eg  — — u"'  in  o i - o«  o*  — co  i _ cc«  — — eg  co  in 

CD  O'  CT*  CO  O'  — eg  eg  CD  CD  O'  CTi  O'  O'  O*  CD  O'  CD  —•  er,  © 

O'  C»  © O'  O O'  CD  — CD  — O'  — — — O'  — 1 CD  — CD  — 

CD  O'  CD  O'  Q O'  O O'  CD  O O'  © O'  O O'  C*  O'  G>  O'  CD 

I II  I I I I I I I 1 1 III 

cr.  — © eg  i'-  Oi  O'  O'  c«  o«  eg  o»  c*  o*  co  oi'  — cd  »£• © 

— in  cc*  rw  is.  v — © O'  ©•  in  o<  o<  c«  ix»  ix  • cd  © r _ in 

in  eg  — eg  — — • rg  c c>  o«  c«  © © © a,  a,  eg  © m •£' 

0 O'  O'  © © O'  O'  O'  O'  O'  in  © o<  o o*  — — u~«  eg  © 

O'  O'  O'  O'  © ©i  O'  O'  O'  O'  O'  O'  O'  O'  O'  O'  O'  © O' 

O'  C*  O'  O'  © O'  O'  © © O'  O'  © O'  O'  O'  © © O'  c<  o» 

++++++++++++++++++++ 

in  ;n  © co  in  ^ in  — eg  — © eg  co  i - in  r . © in  cn  a- 

in  1 - '7i  r-  rv  ij,  — « ^ — . mi  eg  © m oi  ct.  — <£>  in  iv 

cl  eg  — in  fv  i _ cci  in  co  r-  co  'ji'  co  a»  it»  eg  cri  co  — r - cc' 

cd  eg  eg  — © rg  — • eg  — — — — — — — m © lti  — — in 

G»  O'  © O'  C»  © © © O'  © O'  © O'  © O'  C»  O*  O'  © O' 

© © © O'  O'  © © © CD  © O'  O'  O'  O'  © O'  © O'  O'  © 

1 I I I I I 1 I I I I I I I I I I I I I 

rg  r ~ *x'  s7  co  in  *\T  — — eg  eg  eg  ©i  in  © m — in  eg  ©• 

co  zr  ud  in  zr  co  o in  m in  '£'  co  1- 1 - in  co  eg  co  ^ in 

in  co  O'  ud  r v zr  ct.  in  in  in  r - — ir*  ud  in  tt  eg  in  'jd  rg 

eg  ©.  — — ©i  — ©1  ©1  ©1  ©1  — LD  — — eg  f-'.  eg  |s  — '£< 

© ©'  © O'  C>  © © © © O © ©»  © © ©*  © © © ©'  '->» 

C ©•  © © © © © © © © © © ©'  ©'  © © © © © ©' 

+ + + + + + + + + + + + + + + + + + + + 

CT  1 D7  eg  r _ •£.  f.  CTi  'a'  a.  CO  'Ji  UD  K'  in  lv  'X'  © CO  I - 'jJ 

eg  — ■ CD  if»  rg  ^7  © — co  eg  zr  — f - © — I - '£•  e j »_» 

©•  co  ud  ud  — iv  a.  co  1 - ©•  «x<  rn  rg  **T  eg  ud  cr , r--.  eg  co 

— ©i  eg  — rg  iv  od  rg  ^ a>  ud  ? : ^ in  co  — — rg  — ^7 

© ©1  © O'  © O'  © O'  © O'  O'  O'  O'  © O'  © © O'  c«  © 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1111 

©1  Iv  •£•  O',  'jc.  ud  CO  eg  in  — — O'  '£'  UD  in  U"'  UT'  eg  CO  — 

— C'  O'  O'  O'  O'  O'  eg  rg  rg  © — f - © '£•  UD  DT  m eg  '7. 

©i  ©1  O'  O'  O'  O'  O'  O'  O'  O'  eg  © O'  O'  O'  © ' m — eg 

C>  © ©'  O'  O'  O O'  O'  O'  O'  O'  ©'  O'  O'  O’  O'  O'  O'  O'  O' 


+ f + + + + + + + + + + + + + + + + + + 

©•  UD  1 j T Ifi  7 I - co  O'  CD  ^7  UD  CD  «Ji  ^7  '£«  ' D • 1 J • fJ^ 

i - eg  a*  cc  <£.  — a.  u-'  'X*  «'•  j — ud  c*  — ud  'X  D . '£• 

CD  a,  I - u"  eg  IT,  f _ 'X  Z"J.  — T *1  CT'  CD  O'  O'  in  '_J 

Ci-  © C*  O'  O'  — C O'  O'  — Cg  rg  f «•!—■’«  • O'  i—.1 

©«  G © © o ©•  ©’  ©»  ©»  © O'  O*  © ©‘  O'  O'  O'  O'  O'  O' 

I I I I I I I I I I I I I I I I I I I l 

— • rj  CD  eg  co  m 'X*  CC'  I - O'  O'  'X'  O'  CD  eg  u~  CO  DT  CC*  — 
eg  C J a.  © ud  cc  rn  I - eg  UD  in  m O'  I - UD  rn  eg  'X'  I”  ' '£' 
lx  eg  rn  rn  eg  h'  i - in  n m m CO  rvj  C- r - — O'  CO  — fn 
©•  O'  O'  O'  © O'  O'  O'  O'  O'  rn  © — — — — eg  >^_j  u_< 
LT:  O O'  O © O © © O © O O © © © © © ' — 1 

IxJ  J j J • j J • • J • ^ _•  | J|  _•  ( _• 

^ + + + 4-  + + + + “ + + :t::?+:;::  + -t-  + :t: 

h-  t eg  'x>  rn  eg  r _ it,  u-'  — — co  r - rn  *T  © © '£•  <Ti  u"  o, 

<r  tj-  ud  o-  'T  eg  t — co  '£'  ud  co  :o  fs.  <t.  eg  — U”»  m 'J,  IT, 

> «x  '£•  u~  o f J Oi  '£•  r*“'  — u-'  — ' O'  r - r-  *J|  eg  'Ji  CO  — '£•  Di 

• j — i w I w C'j  I - CC'  ee.  eg  i.t,  cc  u~  r - i - CO  u“'  I - U-'  1 ~ t , — — 

Q_  LlI  O'  O'  O'  O'  © O'  — O'  O'  © O'  O'  O'  O'  O'  O'  O'  O'  — — 

LlJ  * * £ 

f i o © O'  O'  © O'  O'  O’  O'  O'  O'  O'  O'  O'  O'  O'  '_J  '— ' ' 

Z © — eg  r*-»  rr  ud  •£•  I - CO  © — eg  rn  T ud  •£•  © © 

cl 


TABLE  10.  (concluded) 
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TMJLF  11.  (continupd) . 


TABLE  11.  (concluded) 


Runs  8 through  11  in  Table  10  are  estimates  from  rudder 
doublet  input  tine  histories  for  Flight  1. . Again  standard 
deviations  are  reduced  for  all  stability  derivatives  and  the 
quality  of  rudder  derivative  estimates  is  much  improved  as 
expected.  Table  11  shows  that,  as  with  aileron  inputs, 
variances  are  reduced  for  soae  derivatives  when  separated 
according  to  initial  rudder  deflection  (2  positive  , 2 
negative) . This  evidence  of  nonlinear  response  suggests 
that  further  research  of  the  analytical  nodel  is  needed. 

Runs  12  through  15  in  Table  10  are  estimates  froa 
side-force  surface  .input  time  histories  for  Flight  1. 
Because  of  a side-force-surf ace-to- aileron  interconnect  and 
a side-f orce -surface- to-rudder  interconnect  used  on  this 
flight  to  better  simulate  pure  lateral  translations, 
side-force  surface  control  derivatives  were  estiaated 
holding  the  aileron  and  rudder  derivatives  fixed  at 
previously  estimated  values. 

Runs  16  through  20  in  Table  10  are  estimates  froa 
side-force  surface  input  time  histories  for  Flight  2.  No 
control  interconnects  were  used  on  this  flight.  The  single 
positive  and  two  negative  side-force  surface  pulses  (15,  19, 
and  20)  produced  larger  state  perturbations  (and  thus  higher 
signal  to  noise  ratios)  than  did  the  single  positive  and  two 
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negative  doublets  (16,  17,  and  18).  Cramfer-Rao  bounds  for 

both  types  of  inputs  are  similar  however.  Table  11 

separates  estimates  by  direction  of  initial  input  and  again 

the  variance  in  some  derivatives  is  reduced. 

Chapter  h discusses  the  guestionable  C derivative 

3 n 

6a 

estimate.  The  estimate  is  5 times  as  large  as  the  analog 

matching  result  and  is  opposite  in  sign.  Table  12  lists  the 

estimates  from  two  aileron  pulse  time  histories  where 

aileron  control  derivatives  were  fixed  at  analog  matching 

values.  In  comparing  the  estimates  with  those  for  runs  1 

and  0 from  Table  10,  C , C , C , and  C estimates 

1r  np  nr 

doubled  in  magnitude  although  C did  decrease  towards  the 

°P 

analog  matching  result. 
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TABLE  12. DERIVATIVE  ESTIMATES  FOR  RUNS  HOLDING  C FIXED  AT 
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APPENDIX  B 

TABULATED  LONGITPDINAL  DERIVATIVE  ESTIMATES 

Longitudinal  stability  and  control  derivatives  were 

identified  by  the  Modified  Maximum  Likelihood  Estimation 

program  using  elevator  doublet  input  time  histories  from 

Flight  2.  Table  13  lists  the  estimates  and  lover  bounds  to 

the  standard  deviations  (Cramer-Rao  bounds)  for  initial  runs 

using  4 negative  initial  input  doublets  (runs  1 through  4) 

and  3 positive  initial  input  doublets  (runs  5 through  7). 

Listed  also  are  the  averages  and  standard  deviations  for  all 

7 runs,  the  negative  runs  and  the  positive  runs. 

Table  14  lists  the  estimates  from  2 time  histories  (4 

and  7)  for  varying  angle  of  attack  upvash  correction 

factors.  Angle  of  attack  measurements  were  divided  by 

correction  factors  of  1,  1.5,  and  1.9.  C appears  to  be 

\ 

the  most  sensitive  of  all  estimated  derivatives  to  errors  in 

the  upvash  correction.  The  variance  in  C is  not  reduced 

Zq 

however  by  changing  upvash  correction  factor. 
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CZ  DER 1 VAT IVES 

RUN  ALPHA  0 

DE 

1-3.28187+0.15433  20.1583+ 

8 . 956,-0 . 82853+0 . 08952 

2-3.88989+0.14645  4.1088+ 

6.858’  -1.1 1753+0. 10091 

3-4. 07454+0. 15917  2.2472+ 

6.158  -1 .6839940. 14265 

4-3.80012+0.12230  9.2276+ 

4.883  -1.1 1669+0.68668 

5-3.86492+0.13668  19.7470+ 

6.482  -1.4O569+0.09275 

6-3.28039+0.13495  54.5444+ 

8.528  -0.66660+0.O9959 

7-3 . 62405+0 . 08399  22 . 5053+ 

4.218  -1.22202+0.86589 

MEAN  AND  STD  DEV  BASED  OH 

7 SAMPLES 

7-3 . 65654+0 . 27307  1 8 . 9338+ 

16.368  -1 . 14958+0.31653 

7.5* 

86.57  27.5* 

MEAN  AND  STD  DEV  BASED  ON- 

4 SAMPLES 

4-3 . 70660+0 . 28542  9 . 9350+ 

6.965  -1.18818+0.31275 

7.7'/. 

77 . 9*  26 • 3* 

MEAN  AND  STD  DEV  BASED  OH  + 

3 SAMPLES 

3-3. 58979+0. 239S4  32.2656+ 

15.734  -1.09810+0.31420 

6.7* 

48 . 9*  28 . 6* 

CM  DERIVATIVES 

RUN  ALPHA  0 

DE 

1 -0 . 38583+0 . 02432- 1 8 . 2562+ 

0.907  -1 .28474+0.04082 

2-0 . 20756+0 . 02O28-26 . 7888+ 

1.446  -1.56401+0.06154 

3-0 . 35 77 4+0 . 02085-26 . 3030+ 

1.731  -1.83234+0.08791 

4-0 . 46323+0 .01841 -20 . 9779+ 

1.063  -1.40250+0.04626 

5-0.52380+0.021 10-10.2095+ 

0.607  -0.86886+0.02361 

6-0.51788+0.02297-1 1 .4985+ 

0.674  -0.91529+0.02880 

7-0.39200+0.01098- 16. 1359+ 

0.566  -1.07576+0.02407 

MEAN  AND  STD  DEV  BASED  ON 

7 SAMPLES 

7-0 . 40675+0 . 10127-18. 5C|56+ 

6.083  -1 .27764+0.32613 

24.9* 

32.7*  25.5* 

MEAN  AND  STD  DEV  BASED  ON- 

4 SAMPLES 

4-0 . 35359+0 . 09274-23 . 08 1 3+ 

3.600  -1.52090+0.20533 

26 . 2* 

15-6*  13.5* 

MEAN  AND  STD  DEV  BASED  ON + 

3 SAMPLES 

3-0 . 47763+0 . 06059- 12.61 46+ 

2.545  -0.95330+0.08864 

12.7* 

20-2*  9.3* 

TABLr,  13.  LONGITUDINAL  DERIVATIVE  ESTIMATES  AND 

CR  A M £p - R AO  BOUNDS. 
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TABLE 


CZ  DERIVATIVES 

FACTOR  ALPHA  0 

1 tj  - 3 . 68? 5 1 +6 . 080 11  44 . 3280+ 
-3. 4? 139+0. 09503  96.2676+ 

DE 

5.043  -2.45525+0.07534 
7.589  -1.86384+0.08109 

1 . 5 -4.68977+0.09312-19.8944+ 

-4 . 62499+0 . 1 1 889  9 . 7629+ 

3.566  -2 . 27630+0 .066861 
5.2O0  -1 .86261+0. 07459' % 

1 . 9 -5.26684+0. 1 1 593-43 . 7695+ 

-5.31 585+0 . 1 433 1 - 20 . 0994+ 

3.496  -2.24196+0.07280 

4. 726  - 1 . 89297 +0. 07756 

CM  DERIVATIVES 

FACTOR  ALPHA  0 

1.0  -0 . 36369+0 . 0 1 230- 1 2 . 2338+ 

-0 . 343 1 2+0 . 0 1 404- 1 0 . 8432+ 

DE 

0.516  -0.9C072+0.00609 
0.756  -0.81520+0.00924 

! c -0.50713+0.01797-10.8547+ 

’ 3 -0.47559+0.02214  -9.1165+ 

0.486  -0.91O96+0.00633 
0.743  -0.82046+0.01019 

i q -0.61282+0.02505  -9.9973+ 
-0 . 5736 1 +0 . 02987  -8 . 1 40 1 + 

0 . 546  -0.91 369+0 . 00.753 

0 . 80 1 -0.82043+0.011 39 

14.  LONGITUDINAL  DERIVATIVE  ESTIMATES  FOR  VARYING 
UFWASH  CORRECTION  FACTORS  (1.0,  1.5,  AND  1.9). 
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APPENDIX  C 


A 

a 

y 

a 

n 

B 

C 

1 

C 

m 

C 

n 

C 

X 

C 

Y 

C 

Z 

c 

c 

o 
D 1 

D2 

f 


G,H 


9 

I 

X 


I 

Y 


NOMENCLATURE 

fundamental  matrix 

lateral  acceleration,  positive  right  (g) 
normal  acceleration  , positive  down  (g) 
control  effect  matrix 

dimensionless  rolling  moment  coefficient 
dimensionless  pitching  moment  coefficient 
dimensionless  yawing  moment  coefficient 
dimensionless  axial  force  coefficient 
dimensionless  side  force  coefficient 
dimensionless  normal  force  coefficient 
vector  of  parameters  to  be  estimated 
S priori  value  of  c 
instrument  error  weighting  matrix 
a priori  parameter  estimate  weighting  matrix 
vector  of  nonlinear  equations  of  motion 
observation  matrices 

2 

acceleration  due  to  gravity  (ft/sec  ) 

moment  of  inertia  about  the  longitudinal  axis 
2 

(slug-ft  ) 

2 

moment  of  inertia  about  the  lateral  axis  (slug-ft  ) 
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moment  of  inertia  about  the  normal  axis  (slug-ft  ) 


cross-product  of  inertia  about  the  longitudinal 

2 

and  lateral  axes  (slug-ft  ) 


cross-product  of  inertia  about  the  longitudinal 

2 

and  normal  axes  (slug-ft  ) 


cross-product  of  inertia  about  the  lateral 

2 

and  normal  axes  (slug-ft  ) 


cost  function 


overall  a priori  weighting  factor 


livelihood  function 


rolling  moment  divided  by  moment  of  inertia  about 

2 

the  longitudinal  axis  (rad/sec  ) 


pitching  moment  divided  by  moment  of  inertia  about 

2 

the  lateral  axis  (rad/sec  ) 


aircraft  mass  (slugs) 


yawing  aoient  divided  by  aoment  of  inertia  about 

2 

the  normal  axis  (rad/sec  ) or 


number  of  sampling  instants 


measurement  errors  vector 


roll  rate,  perturbation 


pitch  rate,  perturbation 


acceleration  transformation  matrix 


C - 2 


i 


r 

5 
u 

u 

V 

V 
K 

V 

w 

X 

X 

X 

V 

y 

i 

z 

7 

Z 

a 

6 

6a 

8e 

6r 


yaw  rate,  perturbation 

2 

reference  area  (ft  ) 
longitudinal  velocity,  perturbation 
control  deflections  vector 
velocity  (ft/sec) 
lateral  velocity,  perturbation 
aircraft  weight,  (lb) 
normal  velocity,  perturbation 
disturbances  vector 

2 

longitudinal  force  divided  by  lass  (ft/sec  ) 
longitudinal  position 
states  vector 

side  force  divided  by  lass  and  velocity  (rad/sec) 
lateral  position 
output  vector 

Donal  force  divided  by  mass  and  velocity  (rad/sec) 

normal  position 

observation  vector 

angle- of-attack 

sideslip  angle 

aileron  deflection,  (right  minus  left)/2,  (deg) 

elevator  deflection  from  trim  (deg) 

rudder  deflection  from  fuselage  centerline  (deg) 
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esf  sideforce  panel  deflection,  (right  plus  left)/2, 

e pitch  angle  (deg) 

i roll  angle  (deg) 

4>  yaw  angle  (deg) 

V gradient 


(deg) 
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